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Abstract. The main goal of this paper is to demonstrate the usefulness of 
certain ideas from System Theory in the study of problems from complex 
analysis. With this paper, we also aim to encourage analysts, who might not be 
familiar with System Theory, colligations or operator models to take a closer 
look at these topics. For this reason, we present a short introduction to the 
necessary background. The method of system realizations of analytic functions 
often provides new insights into and interpretations of results relating to the 
objects under consideration. In this paper we will use a well-studied topic from 
classical analysis as an example. More precisely, we will look at the classical 
Schur algorithm from the perspective of System Theory. We will confine our 
considerations to rational inner functions. This will allow us to avoid questions 
involving limits and will enable us to concentrate on the algebraic aspects of 
the problem at hand. Given a non-negative integer n, we describe all system 
realizations of a given rational inner function of degree n in terms of an 
appropriately constructed equivalence relation in the set of all unitary (n+l) x 
(n + l)-matrices. The concept of Redheffer coupling of colligations gives us the 
possibility to choose a particular representative from each equivalence class. 
The Schur algorithm for a rational inner function is, consequently, described 
in terms of the state space representation. 
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NOTATION: 

T is the unit circle in the complex plane: T = {tGC: = 
D is the unit disc in the complex plane: B = {z G C : |z| < 1} 
is the exterior of the unit circle: = {z : 1 < \z\ < oo}. 



2 



Bernd Fritzsche, Victor Katsnelson and Bernd Kirstein 



dJlpxq is the set of all p x g (p rows, q columns) matrices with complex entries . 
/„ - the identity n x n matrix. 
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0. Introduction 

Up until the 1960s System Theory suggested that a system be considered only 
in terms of its input and output. A system was treated as a 'black box' with 
input and output terminals. Associated with each system was an 'input-output' 
mapping, considered to be of primary importance to the theory at the time. This 
approach, however, did not take the internal state of the system into account. It is 
to be assumed that an input signal will, in some way, influence the internal state 
of a system. Nevertheless, there was little discussion of the relationship between 
input and the inner state of a system until the introduction of State Space Sys- 
tem Theory. This theory not only incorporated input and output spaces, serving, 
respectively, as 'domains' for input and output 'signals', but also a 'state space'. 
This state space was introduced to describe the interior state of the system. 

State Space System Theory (both the linear and general non-linear varia- 
tions of the theory) was developed in the early 1960s. Two names closely asso- 
ciated with the early development of this theory are those of R. Kalman and 
M.S. Livshitz. Kalman's first publications pertaining to State Space System The- 
ory include |KalH [Kal2|, IKal3] . The monograph [KFAJ summarizes these papers, 
among others. R. Kalman's approach to State Space System Theory was from 
the perspective of Control Theory. This approach suggested that the questions of 
a system's controllability and observability be given the most attention. Control 
Theory does not, however, put much emphasis on energy relations and, as a result, 
Kalman's work does not address the subject of energy balance relations (Kalman's 
approach to System Theory was abstract. He develops the theory over arbitrary 
fields, not specifically over the field of complex numbers). In Kalman's theory, one 
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first starts from the input-output behavior (i.e. transfer function) and then con- 
structs the state operator. In Livshitz's theory, the reverse approach is used: The 
characteristic function (which is the analogue of the transfer function) is produced 
from the main operator (which is the analogue of the state operator). Kalman's 
theory is mainly finite-dimensional and afhne, whereas Livshitz's theory is mainly 
infinite-dimensional and metric. It took some decades before the connections be- 
tween these two theories were discovered in the 1970s. Among others, Dewilde 
[Dewlj IDew2) and Helton [Hell IHe2| IHe3| produced much of the work leading to 
this discovery. The connections between the two approaches were made explicit in 
the monograph jBGKj . 

M. S. Livshitz, a pioneer in the theory of non-self-adjoint operators, chose 
to approach State Space Theory from the perspective of Operator Theory. For a 
particular class of non-self-adjoint operators, Livshitz was able to associate each 
operator of this class with an analytic function in the upper half-plane or unit 
disc. These analytic functions were called 'characteristic functions'. Livshitz was, 
furthermore, able to determine a correspondence between the invariant subspaces 
of a linear operator and the factors of its characteristic function (See |Liv3j and 
references within ILiv3| l. Using the framework provided for by these results for 
characteristic functions, Livshitz constructed triangular models of non-self-adjoint 
operators (Triangular models were later partially supplanted by functional models. 
See jSzNFo] ) . Following this, Livshitz focused on questions in both mathematics 
and physics. Oscillation and wave propagation problems in linear isolated systems 
are related to self-adjoint operators. In the mid-1950s M.S. Livshitz began to look 
for a physical example to which his theory of non-self-adjoint operators could be 
applied. This lead him to consider a number of concrete linear systems. These 
systems were not isolated systems, but were such that they allowed for the ex- 
change of energy with the 'external world'. The model of the dynamical behavior 
of a system of this type makes use of an operator and this 'principal' operator 
is, in general, non-self-adjoint. The energy exchange of the system is reflected in 
the non-self-adjointness of the operator. Livshitz worked on problems involving the 
scattering of elementary particles (See |Liv4| . [Liv5j . [BrLi] ) . problems in electrical 
networks (See [LiFl] ) and questions dealing with wave propagation in wave-guides 
(See |Liv6| ) . It was at this juncture that the notion of an 'operator colligation' 
(also common are the terms 'operator node' and 'operator cluster') was introduced 
to provide further clarity. An operator colligation consists of the aforementioned 
'principal' operator, but also 'channel' spaces and 'channel' operators, of which 
the latter two objects describe the non-self-adjointness of the 'principal' operator. 
The introduction of this concept allowed a characteristic function to be associated 
with an operator colligation, as opposed to its respective 'principal' operator (See 
[BrLi] ■ [Br| . |LiYa| and references therein). At much the same time, the concept of 
an 'open system' was then being established (What Livshitz then referred to as an 
'open system' was, in essence, what is now known as a stationary linear dynamical 
system). Livshitz first introduced the notion of an 'open system' in his influential 
paper, [Liv9| (See Definition 1 on p. 1002 of the original Russian paper [Liv7] or 
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p. ??? of the English translation in the present volume). To each system there is 
an associated colligation and in |Liv7] it is shown that a system's transfer oper- 
ator coincides with the characteristic function of the system's colligation. [Liv8| 
introduces the operation of coupling open systems as well as the concept of clos- 
ing coupling channels. jLivSj furthermore introduces the 'kymological resolution' 
of an open system, i.e. the resolution of this system into a chain of simpler cou- 
pled open systems. These simpler systems correspond to the invariant subspaces 
of the 'inner-state' operator of the original open system. To emphasize that the 
notion of an open system is closely related to oscillations and to wave-propagation 
processes, Livshitz uses the terminology 'kymological', 'kymmer' and 'kymmery', 
derived from the Greek word ^KVfia\ meaning 'wave'. We quote from page 15 of 
the English translation of [Liv9' and mention that: "the appropriate representation 
of an open system, transforming a known input into a known output, depends on 
which are known and unknown variables, so that the concept of an open system 
is 'physico-logical' rather than purely physical in nature." 

The relevant theory of open systems and operator colligations, as developed 
by Livshitz and other mathematicians, is presented in the monographs |Liv9| . 
[LiYa] and [Brj . Chapter 2 of the monograph |Liv9] deals with the details of the 
kymological resolution of open systems (a concept of which much use is made in 
the following). A detailed presentation of Scattering Theory for linear stationary 
dynamical systems (with an emphasis on applications to the Wave Equation in 
M") can be found in [LaPhij . 

General State Space System Theory, as developed by R. Kalman and M. S. 
Livshitz provides us with the proper setting and the necessary language for the 
further study of physical systems and various aspects of Control Theory. Despite 
the fact that State Space System Theory does not immediately lead to a solution 
of the initial physical or control problem, it does lead to some interesting related 
questions (mostly analytic). It should, furthermore, be noted that general State 
Space Theory's importance extends beyond its significance within Control Theory 
and when applied to physical systems. M. S. Livshitz was very likely the first to 
understand that this theory had wide-reaching applications within mathematics, 
e.g. in Complex Analysis. 

Analytic functions can be represented or specified in many ways, e.g. as 
Taylor-series, by decomposition into continuous fractions, or via representations as 
Cauchy or Fourier integrals. In the early half of the 1970s an additional method for 
representing an analytic function was introduced, namely the method of 'system 
realization'. This theory has its origins in Synthesis Theory for linear electrical net- 
works, the theory of linear control systems and the theory of operator colligations 
(and associated characteristic functions). M. S. Livshitz established the Theory of 
System Realizations and L. A. Sakhnovich, a former Ph.D. student of Livshitz's, 
later made further important progress in the theory (See |Sakhlj and also [Sakh2j 
for a more detailed presentation of these results). L. A. Sakhnovich studied the 
spectral factorization of a given rational matrix-function R, where both R and the 
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inverse function are transfer functions corresponding to Unear systems (oper- 
ator coUigations). Unfortunately, the paper |Sakhl) did not garner the attention it 
deserved at the time. L. A Sakhnovich's factorization theorem is a predecessor to 
a fundamental result due to Bart/Gohberg/Kaashoek/van Dooren [BGKVj . which 
was remembered as Theorem 2 in the Editorial Introduction to [CWHF] . where 
one can also find a detailed account of the history of the state space factorization 
theorem. 

Our goal is not to provide a comprehensive survey of the history of System 
Theory, so that we have focused on the period leading up to the mid-1970s (with 
particular emphasis on the contributions of M. S. Livshitz and his co-workers). 
His work on open systems was unkown in the western world until his monograph 
[Liv9| was translated in 1973. His fundamental papers |Liv7j and [Liv8j remained 
untranslated up until this memorial volume. 

The subsequent development of the Theory of System Realizations is gener- 
ally associated with the name I. Gohberg, who produced and inspired much in the 
way of new work and results for this theory and its applications. As a result, the 
theory experienced a period of accelerated growth, beginning in the late 1970s. 
Published in 1979, the monograph [B_GkJi| dealt with general factorizations of a 
rational matrix-functions as well as with the Wiener-Hopf factorization of rational 
matrix-function, where, in both cases, this function is a transfer function for a 
linear system (operator colligation). 

I. Gohberg and his co-workers have shown that State Space Theory has a 
much wider range and goes far beyond System Theory and the theory of operator 
colligations. We list a few topics to which State Space Theory can be applied: 

1. Methods of factorization of matrix- and operator- valued functions; solutions 
of Wiencr-Hopf and singular integral equations. 

2. Interpolation in the complex plane and generalizations. 

3. Limit formulas of Akhiczcr/Kac/Widom type. 

4. Projection methods, Bezoutiants, resultants. 

5. Inverse problems. 

The monograph [BGR| offers a detailed discussion of interpolation problems 
and many other questions. Matrix-function factorization is a tool applied in dis- 
cussions of many other problems as well, e.g. in the theory of inverse problems 
for differential equations and also in prediction theory for stationary stochastic 
processes. If a matrix-function is rational, then this factorization can be attained 
using system realizations. These system realizations, in turn, play a certain role 
in the solution of the original problem (See, for example, [AGj ). The Theory of 
Isoprincipal Deformations of Rational Matrix-Functions (which is, in particular, 
a useful tool for investigating rational solutions of Schlesinger systems) is formu- 
lated in terms of the Theory of System Realizations (See [KaVol| and |KaVo2] . 
For our purposes, the theory developed in |Ka| is most relevant). The current state 
of System Theory, as a branch of pure mathematics, is presented in [Nikj . 



^ N.b. There is now an extended version of this monograph. See |BGKR| . 
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In the present paper we show how the Schur algorithm for contractive holo- 
morphic functions in the unit disc can be described in terms of system reahzations. 
In the foUowing, we consider only rational inner functions, which allows us to avoid 
questions involving limits and enables us to concentrate on the algebraic aspects 
of the problem at hand. At first glance the formulas here presented might seem 
rather complicated and, to some degree, less than intuitive. This is, however, from 
the perspective of System Theory, not the case. The aforementioned formulas serve 
as the function-theoretical counterpart to Livshitz's kymological resolution as ap- 
plied to the system (represented by the original inner function) corresponding to 
the cascade coupling, i.e. the Redheffer coupling, of open systems. The elementary 
open systems, which make up this cascade (or chain) correspond to the steps of 
the Schur algorithm. 

This paper is organized as follows. In Section[Tl we state some facts relating to 
rational inner functions. In Section[2l we discuss some aspects of the classical Schur 
algorithm. Section [3] is devoted to a short introduction to operator colligations 
and their characteristic functions, where particular attention is paid to finite- 
dimensional unitary colligations. The characteristic functions of finite-dimensional 
unitary colligations are shown to be rational inner matrix-functions (See Theorem 
13. Sp . Theorem 13.61 shows that an arbitrary rational inner matrix- function can, 
on the other hand, be realized as a characteristic function of a finite-dimensional 
minimal unitary colligation. The scalar rational inner functions of degree n are just 
the finite Blaschkc products of n elementary Blaschke factors. The essential facts 
on the realization of scalar inner rational functions of degree n as characteristic 
functions of minimal unitary colligations are summarized in Theorem l3.10l These 
minimal unitary colligations can be equivalently described by equivalence classes 
of minimal unitary (n -f- 1) x (n + l)-matrices. A proof for Theorem 13.101 can be 
found in the Appendix at the end of the paper. 

The main objective of this paper can be described as follows. The application 
of the Schur algorithm to a given rational inner function s(z) of degree n produces a 
sequence Sfc(z) , fc = 0, 1, . . . , n of rational inner functions with sq{z) ~ s{z) and 
degSi{z) = n — k.ln particular, the function s„(z) is constant with unimodular 
value. In Section [3l it will be shown that each of the functions Sk{z) admits a 
system representation 



in terms of the blocks of some minimal unitary matrix Uk £ 9K(„+i)x(n+i), 



We assume that Uq is given. The goal is to recursively produce the sequence 
matrices J7fe. In other words, the steps of the Schur algorithm have to be described 
in terms of the state space representation. Since the unitary matrices Uk are defined 
only up to an equivalence relation, we have to find corresponding operations for 
the arithmetic of these equivalence classes. 



Sk{z) = Ak + zBk{I ~ zDk) ^ Ck 
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In Section |4] we discuss the means by which the hnear-fractional transfor- 
mation associated with the Schur algorithm can be described in terms of the 
input-output mapping of linear systems. The Redheffer coupling of linear systems 
will be introduced as a useful tool in these considerations. 

In Section [5l the Redheffer coupling of linear systems will be translated into 
the language of unitary systems. 

In Section[Sl we apply the concept of Redheffer couplings of colligations to the 
linear-fractional transformation associated with the inverse of the Schur algorithm. 
In so doing, we will describe the 'degrees of freedom' of unitary equivalence. A 
closer look shows us that amongst all the unitary matrices which realize the desired 
system realization, there are some distinguished by the fact that they are, in a 
sense, associated with the concept of Redheffer coupling. 

In Section [71 the basic step of the Schur algorithm will be described in the 
language of colligations. This requires that we solve a particular equation for uni- 
tary matrices, suggested by the results of Section [51 The solution to this matrix 
equation is given in Theorem l7.1l Together with Lemma [7.2l Theorem l7.2l describes 
the basic step of the Schur algorithm in terms of system representations. 

The investigations of Section [51 show that a certain normalization procedure 
has to be performed at every step of the Schur algorithm if the Schur algorithm is 
to be dealt with in the language of system realizations. We consider the degrees of 
freedom for this normalization procedure. It turns out that we can use these degrees 
of freedom to make the normalization procedure a one-time-procedure, so that 
it might be dealt with during preprocessing for further step-by-step recurrence. 
A one-time-normalization of this kind is related to the reduction of the 'initial' 
colligation matrix to the lower Hessenberg matrix. 

In Section [31 we will be well-positioned to present the Schur algorithm in 
terms of unitary colligations representing the appropriate functions. 

In Section [TUl we express the colligation matrix in terms of the Schur param- 
eters. 

In the final section (Section [TT|) we discuss some connections between the 
present work and other work relating to the Schur algorithm as expressed in 
terms of system realizations. In particular, we discuss the results presented in 
Alpay/Azizov/Dijksma/Langer [AADL] and Killip/Nenciu [KiNej . 

1. Rational Inner Functions 

We say that a function s : D ^ C, where s is holomorphic in D, is 
contractive if 

\s{z)\ < 1 for every ^ G D. 
A contractive function s is called an inner function if 

\s{t)\ = 1 for every t e T. 
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In the following we consider rational inner functions, so that s{t) is defined for 
every t <E T. 

A rational function is representable as a quotient of irreducible polynomials 
and we call the order of the highest-degree polynomial the degree of the rational 
function. 

If a rational function s is an inner function, then the degree of its numerator 
and the degree of its denominator are equal. 

An inner rational function s is representable as a finite Blaschke product, i.e. 
in the form 

n (") 

l<fe<n " 

Zi, . . . , Zn are points in D, or, in other words, complex numbers satisfying the 
condition 

|zi| < 1, ... , |z„| < 1, (1.2) 
c is a unimodular complex number, i.e. 

|c| = 1. (1.3) 

Conversely, given complex numbers zi , . . . , z„ and c satisfying the conditions (|1.2p 
and (|1.3p . respectively, the function s in p.ll) is an inner rational function of degree 
n. 

The number c and the set {zi, ... , z„} are uniquely defined by the inner 
function s (the sequence of numbers (zi, . . . , z„)) up to permutation. 

The notions of contractive and inner functions can also be defined for matrix- 
functions: 

We say that a matrix-function S* : B ^ 97lpxp , where S is holomorphic in B, 
is contractive if 

Ip - S*{z)S{z) > for every z G D . 
A contractive matrix-function 5' : B ^ ^pxp , is called an inner function 1:^1 
Ip - S*{t)S{t) = for almost every i e T . 



2. The Schur Algorithm 

In this section, we present a short introduction to the classical Shur algorithm, 
which orginated in Issai Schur's renowned paper, [Schj. In so doing, we will mainly 
emphasize those aspects of the Schur algorithm, which are essential for this paper. 
For comprehensive treatments of the Schur algorithm and its matricial general- 
izations, we refer the reader to jBFKl] . |BFK2| . |Con2| . [PFK] . |S:Meth| and the 
references therein. 

Let s{z) be a contractive holomorphic function in B and 

so = s{0). (2.1) 

^For a contractive holomorphic function S in D, the boundary values S(t) lim S{rt) exist 

r— >i— 

for almost every t G T (with respect to the Lebesgue measure). 
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Then |so| < 1, where |so| — 1 only if s(z) — sq- If |so| < 1, then the function 

= (2.2) 

is well-defined. Moreover, it is contractive holomorphic in D. The function s(z) 
can be expressed in terms of these w(z) and sq" 

s{z) = IR±I^ . (2.3) 

If the function 5(2:) is an inner function, then ioi^z) is also an inner function. If s(z) 
is an inner rational function of degree n, then ijo{z) is an inner rational function of 
degree n — 1. 

Conversely, if 07(2:) is an arbitrary contractive holomorphic function in D 
and So is an arbitrary complex number satisfying the condition |so| < 1, then 
the expression on the right-hand side of (|2.3p defines the function s(z), which is 
holomorphic and contractive in D. Furthermore, if uj(z) is an inner function, then 
s{z) is an inner function as well. 

DEFINITION 2.1. 

/. We call the transformation s{z) 1 — > '^{z), defined by (12. 2p . where 
sq = s(0), the (direct) Schur transformation. 

//. We call the transformation uj{z) 1 — > s{z), defined by (|2.3p . where Sq is a 
given complex number, the inverse Schur transformation. 

The correspondence s{z) (s(0), i-^{z)) describes the elementary step of 
the Schur algorithm. 

The Schur algorithm is applied to a holomorphic function s(z), which is 
contractive in D. This algorithm inductively produces the sequence (finite or in- 
finite) of contractive holomorphic functions Sfe(z) in D and contractive numbers 
Sfe = Sfe(O), fc = 0, 1, 2 , . . . . The algorithm terminates only if s(z) is a rational 
inner function. Starting from s(z), we define 

sq{z) = s{z), So = so(0) . 

If the functions Si{z), i = 0, 1, . . . , k are already constructed and 
|sfe(0)| < 1, then we construct the function Sk+iiz) as follows: 

f \ 1 Sfe - Sk{z) 

Sfe+i(z) = -- — — , Sfe+1 = Sfc+i(0) . (2.4) 

Z 1 - Sk[Z)Sk 

If s{z) is not a rational inner function, then the algorithm does not terminate: 
On the k-th step we obtain the function Sfe(z), for which 
Sfc(0)| < 1, so that we can construct the function Sk+i{z) and still have 
|sfe+i(0)| < 1. 

If s(z) is a rational inner function of degree n, then we can define the functions 
Si{z) for i = 0, 1, . . . , n such that 

deg Si{z) — n — i, z = 0, 1, . . . , n . 
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The numbers Si = Si(0) satisfy the conditions 

\si\ < 1, i = 0, 1 . . . , n - 1 . 

However, in this case 

\^n\ — I5 ^n(^) = ^71- 

So, for k = n the numerator and the denominator of the expression on the right- 
hand side of (|2.4p vanish identically. The function Sn+i{z) is thus not defined and 
the Schur algorithm terminates. 

The numbers Sk = Sfc(O) are called the Schur parameters of the function s{z). 

If s(z) is not an inner rational function, then the sequence of its Schur parame- 
ters is infinite and these parameters Sk satisfy the inequality 
\sk\ < 1 for all A; : < A; < 00. If s{z) is an inner rational function 
with degs(z) = n, then its Schur parameters Sk are defined only for 
k — 0, 1, ... ,n and 

|sfe|<l, A: = 0, 1, ... , 71-1, |s„| = l. (2.5) 

Conversely, given complex numbers sq, si, ... , s„ satisfying the conditions (|2.5p . 
one can construct the inner rational function of degree n, having Schur parameters 
So, si, . . . , s„. This function s(z) can be constructed inductively: First, we set 

If the functions Si{z) for i = n, n — 1, . . . , k are already constructed, then we set 

Sk-l + ZSkiz) 

Sk-l{z) = 1-T ■ 

1 + ZSk~l Sk{z) 

In the final step we construct the function so{z) and set 

s{z) = so{z). 

Thus, there exists a one-to-one correspondence between rational inner func- 
tions of degree n and sequences of complex numbers {sk}o<k<n satisfying the con- 
ditions (12. 5p . 

3. The System Representation of a Rational Inner Function. 

Contractive holomorphic functions appear in several roles. In particular, such func- 
tions appear in Operator Theory as the characteristic functions of operator colli- 
gations. The notion of an operator colligation is closely related to that of a linear 
stationary dynamical system. There is a correspondence between the theory of 
operator colligations and the theory of linear stationary dynamical systems. The 
concepts and results of one theory can be translated into the language of the other. 
There are interesting connections to be made between these theories. Definitions 
and constructions, which are well-motivated and natural in the framework of one 
theory may look artificial in the framework of the other. In particular, the notion 
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of the characteristic function of a colligation and of the coupling of colligations are 
more transparent in the language of System Theory. 

In this section, the term 'operator' means 'continuous linear operator'. 

DEFINITION 3.1. Let H, 5°"* be Hilhert spaces and U be an operator: 

U : r" ®H 5°"* ©H, (3.1) 



U = 



(3.2) 



Let 

'A B 
C D 

be the block decomposition of the operator U, corresponding to (|3.ip ; 

A : r" -> B-.n^ 5°"*, C : -^n,D: H^H. (3.3) 

The quadruple {£'"'■, H, U) is called an operator colligation. 

and f are, respectively, the input and output spaces of the colligation. 
We call Ti the state space of the colligation and A the exterior operator. We call 
B and C channel operators, while D is referred to as the principal operator of the 
colligation. Finally, we call U the colligation operator. 

If the input and the output spaces and coincide: — = £, we 
call the space £ the exterior space of the colligation and denote the colligation by 
the triple {£, Ti., U) 

DEFINITION 3.2. Let (f™, H, U) be an operator colligation. 
The operator-function 

S{z)^ A + zB{In~ zOy^C (3.4) 

is called the characteristic function of the colligation (£'", f ?i, U). 

The function S{z) is defined for the z 6 C where the operator (I-h — zD)^^ 
exists. The values of S are operators acting from into £°"* . 

REMARK 3.1. The function S{z) is defined and holomorphic in some neighbor- 
hood of the point z = 0. Furthermore, S{0) — A. 

The notion of a colligation's characteristic function draws on the framework 
of the theory of linear stationary dynamical systems (LSDS). (The theory of open 
systems, in the terminology of M.S.Livsic). The theory of LSDS, which we are 
dealing with is not a 'black box theory', where only the input signals, output 
signals and the mapping 'input ^output' are considered. The theory of LSDS also 
takes 'interior states' of the system into account. The input and output signals are 
described (in the discrete time case, where the index fc serves as time) by sequences 
\}Pk\^<k<oa and {V'/c}o<fc<oo of vectors belonging to some Hilbert spaces £™ and 
gowt ^|.]-^g jyjpy^ ^]-^g output spaces of the system). The 'interior states' are 
described by vectors /i of a Hilbert space Ti, called the state space of the system. 
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The dynamics of a linear stationary system is described by the hnear equa- 



tions 







'a 


B 






hk+1 




c 


D 




hk 



A: = 0, 1, 2, 



(3.5) 



where the operators A, B, C, D do not depend on k ('time') and are defined in 

(ESI). 

It is natural to consider the four operators A, B, C, D as blocks of the 'uni- 
fied' operator, say U as in (|3.2p . from the space £^'^®H into the space f °"*©7i. The 



operator colligation (£*", H, U) then corresponds to the LSDS 

Given the sequence {'Pk}o<k<m and the initial value /iq, the system (|3.5p uniquely 
determines the sequences {'0fc}o<fc<m and {/ife}o<fc<m+i- In the case /iq = 0, 



TO > 1 . 



(3.6) 



V'o = Alpq, -(/),„ = Aip„i + 2^ BD 'Ctpm-k-i ■ 

l<fe<m-l 

The relation p.6p can be considered as the description of the evolution of the LSDS 
p.Sp in the time domain. The description of the evolution is, however, especially 
transparent in the frequency domain. Since the considered sequences are unilateral, 
the Fourier transforms of these sequences are the (formal) power series 



ip{z) 



0<fc<oo 



•^kZ 



^(z) 



0<fc<oo 



■tpkz'' , h{z) 



0<fc<oo 



hkz" 



(3.7) 



The complex variable z can be interpreted as the frequency. Under the extra 
assumption that /iq = we can rewrite (|3.5p in terms of the Fourier representations: 







'A B 




ip{z) 


z-^h{z) 




C D 




h{z)_ 



(3.8) 



From 



we obtain 



i){z)^ Atp{z)+Bh{z), (3.9a) 

h{z) = z{I - zD)-^Clp{z) . (3.9b) 

Eliminating h{z), we get 

i;{z) = S{z)^{z) , (3.10) 

where S{z) is expressed in terms of the matrix p.2p as in (|3.4p : 

S{z)^ A + zB{I ^ zD)-^C . 

The operator function S{z) describes the input-output mapping corresponding to 
LSDS (1^ . 

DEFINITION 3.3. In the framework of System Theory, the function S{z) in ((XTU)) 
is called the transfer matrix of the LSDS (|3.5p . 



In the theory of operator colligations the operator function S{z) is called 
the characteristic function, while in the theory of LSDS it is called the transfer 
function. This notion, however, makes more sense in the theory of LSDS. Along 
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with the input-output mapping described by the transfer function S{z), the input- 
state mapping: 

<p(z) h{z) , where h{z) = z {I — zD)^^Cip{z) , 

is also naturally related to the system p.Sp . 

If the dimensions dim £™ and dim f °"* of the input and output spaces are 
finite, then, choosing bases in and we can consider S{z) as a matrix- 

valued function. If, moreover, the dimension dimTi of the state space is finite, 
then S{z) is a rational matrix-function. 

DEFINITION 3.4. The colligation f""*, H, U) is said to be finite-dimensional 

if dinii?'" < oo, dimf°"* < oo and dimTi < oo. 

The dimension dim TC of the state space of the finite-dimensional colligation 
f""*, H, U) is related to the degree of its characteristic function. Here we 
use the notion of the McMillan degree of a rational matrix-valued function as 
it is defined in fMcM' . The notion of the degree of a rational matrix-function is 
discussed in [DuHa and [KaMJ. See also |BGK| . In the case when dim£ = 1, 
i.e. in the case when the considered rational function is scalar (or C- valued), the 
McMillan degree of this function coincides with its 'standard' degree. 

To precisely formulate how the dimension of the state space H. and the degree 
of the characteristic function S{z) are related, we need to introduce the notion of 
a minimal colhgation f""*, H, U). 

DEFINITION 3.5. Let {£"\ f""*, H, U) he a colligation. We define the following 
subspaces of the state space TL: 

H" = clos( Y (D'^-C)^"), 7^° = c1os( Y {D*''B*)£'"'*) , (3.11) 

0<fc<oo 0<fc<oo 

where \/ fk denotes the linear hull of the vectors fk and clos{M) denotes the closure 

k 

of the set M . 

The subspaces Ti.'^ and Ti.° are, respectively, called the controllability and ob- 
servability subspaces of the colligation (£™, T[, U). 

REMARK 3.2. // the state space TL is finite- dimensional, say 
dimTi = n < oo, then it is enough to restrict our considerations in (j3.1ip to 
the linear hull of the vectors {D^C)£"^ and {^Q*^ B*)£°^* with k < n. In this case 
there is no need to make use of the closure in (j3.1ip . 

DEFINITION 3.6. We say that a colligation (f™, f""*, H, U) is controllable if 
and observable ifn° = H. 
We say that a colligation is simple if the sum of the controllability and the 
observability subspaces is dense in the state space, i.e. if 



clos{H'' + H°) 
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We say that a colligation (f™, 5°"*, H, U) is minimal if it is both controllable 
and observable, i.e. if 

^ = 1-1 and re = 1-1. 

THEOREM 3.1. Let f™*, H, U) be a finite- dimensional colligation and let 
S{z) be the characteristic function of this colligation. 

S{z) is then a rational matrix- function, which is holomorphic at z = and 
such that 

deg5<dimH (3.12) 
Equality holds in (j3.12p if and only if the colligation (f™, f""*, 7i, U) is minimal. 

THEOREM 3.2. Let £i and £2 be finite- dimensional spaces and let S{z) be a 
rational function, whose values are operators acting from £1 to £2 and which is 
holomorphic at the point z = 0. 

There then exists a finite- dimensional minimal operator colligation 
(f™, f™*, n, U), dSII]) - ^ - dSS]), with £"' = £1 and f™* = £2, whose charac- 
teristic function Su{z) = A-\- zB{I — zD)^^C coincides with the original function 
S{z). In other words, S can be expressed in the form p.4p . 

DEFINITION 3.7. The representation of a given function S{z) as a characteristic 
function of an operator colligation is called the state space representation of S{z) 
or the state space realization of S{z). If the representative operator colligation is 
minimal, then we say that the state space realization of S{z) is minimal. 

Let us discuss the uniqueness of the state space representation. 



DEFINITION 3.8. Let {£{", £°"\ Hi, Ui) and f^"*, Wa, U2) be two operator 

'A2 B2 

U2 = , (3.13) 

C2 D2 



colligations: 



where 
Ai : £1 



Bi : Hi 



Q : £1 



Hi, Di : Hi 



i = l,2. 



(3.14) 



We consider the colligations (£[", f""*, Hi, Ui) and (f™, £:°"*, H2, U2) to be 
equivalent if invertible operators E™ , E""^* and V: 



t, .t2^t^, h . t2 

exist, such that the intertwining relation 



V:H2=Hi, 







'A2 


B2 




'A, 


Bi 






V 




C2 


D2 




Ci 


Di 




V 



(3.15) 



(3.16) 



holds. 
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Clearly, given two equivalent operator colligations, one of these colligations 
is controllable, observable, simple or minimal if and only if the other colligation 
possesses the same respective property. 

The following result is evident: 

THEOREM 3.3. Let (£f , fj*"*, Hi, Ui) and £2°"*, ^2, U2) be operator col- 
ligations. Assume that these colligations are equivalent, i.e. that the intertwining 
relation l|3.16p holds with some invertible operators i?'", E°"* and V. 

Then the characteristic functions Si{z) and 82(2) of these colligations, 

S^{Z)^A^+ZB^{I-ZD^)-^C^, ^=l, 2, (3.17) 

satisfy the intertwining relation: 

£;°"*S'2(z) = (3.18) 
for all z where Si and S2 are defined. 

Under the extra assumptions that the colligations are minimal and finite- 
dimensional we can show that for Theorem 13.31 the converse assertion also holds. 

THEOREM 3.4. Let 81'^\ Hi, Ui) a«d(£:^", f^*, H2, U2) fee finite-dimensional 
operator colligations. Let Si{z) and S2{z), p.l7p . fee the characteristic functions 
of these colligations. We make the following assumptions: 

1. The functions Si{z) and 6*2 (^) satisfy the intertwining relation p.lSp for all 
z small enough, where E"^ : £3" ~^ ^1" o.'^d E°^* : f™* ~* ^1"* ^''e some 
invertible operators. 

2. The colligations (£1", £^'^\ Hi, Ui) and {£^ , £^'^\ H2, U2) are minimal. 

These colligations are then equivalent, i.e. there exists an invertible operator 
V : 0.2 — > Ti-i such that the intertwining relation p.l6p holds. 

Up to this point, we have not taken advantage of any scalar products that 
may be defined in the input, output and state spaces. From this point forward, 
we will focus more on these scalar products and the benefits they bring when we 
have them at our disposal. In what follows, we consider rational inner functions. 
Operator colligations representing such functions are unitary, finite- dimensional 
operator colligations. 

For convenience, we recall the definition of a unitary operator: 
Let Ci and C2 be Hilbert spaces and T : Ci ^ C2 be an operator. We say 
that T is unitary if it satisfies the following two conditions: 

a) T preserves the scalar product, i.e. 

{Tx, Ty)c2 = {x, y)ci Va; £ £1, y G £1 . 

b) T maps Ci onto C2, i.e. T is invertible. 
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The unitarity property of a linear operator T can also be characterized as 
follows: 

T*T = /£, , TT* = lo, ■ 

DEFINITION 3.9. Let (f™, H, U), be an operator colligation. 

We call {£"^, H, U) a unitary colligation if tfie colligation operator U is a 

unitary operator, i.e. if 

U*U = Is^^^H, UU* = l£o..^^. (3.19) 

DEFINITION 3.10. Let £i and £2 he finite- dimensional Hilhert spaces and let S{z) 
be a rational function whose values are operators acting from £1 to £2 ■ 

The matrix-function S is called an inner function if its values S{z) are con- 
tractive operators for z £ D and unitary operators for t ^ T, i.e. if the conditions 



l£,-S*{z)S{z)>0,l£,-S{z)S*{z)>0, for zeB, (3.20a) 

hi - S*{t)S(t) =0, Is^ - S{t)S*{t) = 0, for teT. (3.20b) 
hold. (In particular, S has no singularities m B U T.) 

REMARK 3.3. Since unitary operators are invertible, £1 -£2 inner functions exist 
only if dim £1 = dim £1 . 

THEOREM 3.5. Let (f™, £°"*, H, U), he a finite- dimensional uni- 

tary colligation and S{z), (j3.4p . be its characteristic function. 
Then the function S{z) is a rational inner function. 

PROOF. The proof of this lemma is based on identity p.Sp . where h{z) is expressed 
in terms of 1^9(2) as in (|3.9bp . Let z and C be such that the operators / — zD and 
/ — C^* are invertible (These operators are invertible if z e D, C e D. Also, since 
the spectrum of the operator is a finite set, the operators / ~ zD and / — <^D 
are invertible for all but finitely many z G T, C G T.) Because the operator U is 
unitary, (|3.8p yields 

(V'(^), V'(C))e- + (<)"' ( h{z), h{0 )h = {viz). viOs + { h{z), h{0 )n , 

or 

(^(z),(^(C))^.„ - {s{zMz),s{CMC)),.., = 

- (1 - zO{{I - zA)-'Cipiz), (I - CA)-iC^(C))^ . (3.21) 

In particular, taking (p{z) = ip' and (p{C) = v" , where (p' and <p" are arbitrary 
vectors in f™, we obtain the equality 



1-Cz 

In the same way we obtain the equality 
Ieo..-S{z)_S*{0 
1-zC 



= C*{I -CD*)-\l - zD)-^C . (3.22) 
uality 

= B{I - zD)-\l -CD*)-^B* . (3.23) 
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Using the identity '^^^-'^^^ '^jj-' '-^^ = {I - (D)-\l - zD)-\ we obtain 

^i^l-^^ ^ B{I ^ (D)-^! - zD)-'C , (3.24) 
s> ^ 

and 

S*{q-S*{z) ^ c*{I -lD*)-^{I -zD*)-^B* , (3.25) 

To get ([5:^ we let C = 2 in (j?:^ - (023): 

-S*{z)S{z) = (1 - \z\^)C*{I ~zA*)-^[I - zA)-^C, (3.26a) 

/£o„t - S'(z)S'*(z) = (1 - Izp) B(/ - zA)-i(/ - zA*)-^B* . (3.26b) 

The inequalities p.20ap follow from equalities p.26p . which hold for all z G D. 
The equalities (I3.26P furthermore hold for all but finitely many z e T. Thus, the 
rational function S{z) is bounded in T, except on a finite set. S therefore has no 
singularities in T and takes unitary values there. 



The following theorem serves as a 'unitary' counterpart to Theorem 13.21 

THEOREM 3.6. Let S{z) be a rational inner function whose values are operators 
acting from £i into £2, where £1 and £2 are finite- dimensional Hilbert spaces. 

Then there exists a finite- dimensional, minimal, unitary operator colligation 
(ff , f™*, n, U), dsn]) - dSJ) - dSSl), with = £1 and f™* = £2, whose charac- 
teristic function Su{z) = A-\- zB{I — zD)^^C coincides with the original function 
S{z). In other words, the function S is representable in the form p.4p . 

DEFINITION 3.11. Let (fj", f^""*, Hi, Ui) and (£^ , f^*, H2, V2) be operator 
colligations, p.l3p - p.l4l) . If these colligations are equivalent (i.e. if they satisfy 
the intertwining relation (I3.16P - (|3.15p ] and each of the operators i?*", V 
is unitary, we say that (f"\ f^*, Hi, Ui) and (f™' '^^2"*' "^2, U2) are unitarily 
equivalent. 

Clearly, if two operator colligations are unitarily equivalent and one of these 
colligations is unitary, then the second colligation is also unitary. 

The following theorem provides us with a 'unitary' counterpart to Theorem 

1X31 

THEOREM 3.7. Let (fj", Hi, Ui) and {£f, £^'^\ W2, J/a) be unitary col- 
ligations, (j3.13p . Furthermore, let these colligations be unitarily equivalent, i.e. 
suppose that the intertwining relation p.l6p holds for some unitary operators 
E'", £;°"* and V. 

The respective characteristic functions Si{z) and 6*2(2) of these colligations, 
p.l7p . then satisfy the intertwining relation (|3.18p with these very same unitary 
operators and . 
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If we, furthermore, assume that both unitary coUigations are simple, we can 
show that the converse to Theorem 13 . 71 also holds. 

The next theorem serves as a 'unitary' counterpart to Theorem l3.4l 

THEOREM 3.8. Let{£l", Hi, Ui) and{Ef, f^"*, H2, Vi) he finite- dimensional 

unitary operator colligations, (j3.13p . Let Si{z) and S2{z), (j3.17p . be the charac- 
teristic functions of , f Hi, Ui) and , ^2"*' "^2, U2), respectively. We 
now make the following assumptions: 

1. The functions Si{z) and S2{z) satisfy the intertwining relation p.lSp for z £ 
ID, where E^" : S"^, E°^*- : £2"* — > ff"* are some unitary operators. 

2. The colligations (fj", f^*. Hi, Ui) and , f""*, H2, U2) are simple. 

The colligations (f™, 1?°"*, Hi, Ui) and {S™ , ^2"*' U2) are then unitar- 
ily equivalent, i.e. there exists a unitarily operator V : H2 — >■ Hi such that the 
intertwining relation (|3.16p holds. 

Let us compare the assumptions of Theorems 13.41 and 13.81 In Theorem 13.41 
we assume that the colligations ff"*, Hi, Ui), i = 1,2, are minimal, however 
it is not assumed that these colligations are unitary. In Theorem 13.81 we assume 
that the colligations are unitary and simple, but we do not explicitly assume that 
these colligations are minimal, because they are, in fact, already minimal. 

THEOREM 3.9. Let 5°"*, H, U) he a finite-dimensional, unitary operator 

colligation. The following statements are then equivalent: 

1. The colligation is simple. 

2. The colligation is minimal. 

3. The colligation is controllable. 

4. The colligation is observable. 

In what follows we deal only with scalar-valued inner functions S{z), i.e. with 
functions whose values are complex numbers. The input space and the output 
space of the unitary colhgation H, U) representing this S{z) can 

be identified with the space C: f ™ = = C. The finite-dimensional state space 
H, with, say AimH = n, can be identified with the space C" (with the standard 
scalar product): H ^ C". With these conventions in place, the orthogonal sums 
© H and £°"* © H can be identified naturally with the space C © C". 
We note that C©C" represents a canonical decomposition of the space C""*"^ 
into an orthogonal sum. We consider the space C"+^ as the set 97l(„_(_i)xi of all (ri+ 
l)-column-vectors, along with the standard linear operations and scalar product: 

(/,.9)=.9*/, /, 5eaH(„+i)xi, (3.27) 

where the asterisk * denotes Hermitian conjugation. 

A unitary operator, U, acting in C"+^ is described by a unitary (1 + n) x 
(1 -|- n)-matrix, which will also be denoted by U. U maps the column-vector / to 
the column-vector Uf, where Uf is the usual matrix product. The decomposition 
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£n+i — c £;n Qf ^j^q space C"+^ suggest that we consider the foUowing block- 
matrix decomposition of U : 

'a B 

C D 



U = 



nxn • 

The matrix entries are considered as operators: 

A:£^£, B-.n-^S, C : £ ^ H , D : H ^ H . 

where 

f-OTixi (=C), ?i = 97l„xi (-C"). 



(3.28a) 
(3.28b) 



(3.29) 



(3.30) 



DEFINITION 3.12. Given a unitary matrix U G 9K(„+i)x(n+i) with block decom- 
position p.28p . we associate the unitary colligation {£, Ti, U) with U. The exterior 
space £ and the state space Ti of this colligation are as in (j3.30p , where the spaces 
C and C" have the standard scalar products. The exterior, principal and channel 
operators A, D, B, C correspond to the block-matrix entries in (|3.28p and satsify 
(13:291) . 

We call this colligation the unitary colligation associated with the unitary matrix 

U. 

Given two unitary colligations associated with unitary matrices U' and U" , 
how do we express that these colligations are unitarily equivalent? The exterior 
spaces of both colligations are 'copies' of the same space C. To identify the exterior 
spaces C of two different colligations, we should specify the unitary operators i?*" 
and for the two copies of C (These operators, E™ and appear in (j3.15p 

and in the intertwining relations ()3.16p and (|3.18p .) We can naturally choose these 
identification operators as the identity operators, i.e. such that each of operators 
i?*" and is represented by the 1 x 1-matrix whose (unique) entry is the 

number 1 (Such operators can be represented by 1 x 1-matrices, where the matrices 
corresponding to E™ and E°"* consist, respectively, of an arbitrary number i/*" 
and v°"* with = 1 and |t/°"*| = 1.) 

With this convention in place, the unitary equivalence of the colligations 
associated with the block-matrices 



U' 



B' 
D' 



{n+i)x{n+i) and U' 



A" 
C" 



B" 
D" 



means that these matrices satisfy the intertwining relation: 



e 9JI(„+l)x(ri+l) 

(3.31) 



"l 


o' 




'A" 


B"' 




'A' 


B'' 




"l 


o' 





V 




C" 


D" 




C 


D' 







V 



(3.32) 



where V € dJlnxn is a unitary matrix. The equality (j3.18p then becomes: 

Si{z) = S2{z). 
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DEFINITION 3.13. We say that the unitary matrices U' E ^(n+i)x(n+i) O'^d 
U" G 9Jl(„+i)x(ri+i); (|3.3ip . are equivalent if there exists a unitary matrix V G 
9KnxTi such that the intertwining relation (j3.32p holds. 



Let U = 



G 9K(„+i)x(n+i) • We now consider the following matrices 



A B 
C D 

associated with the unitary matrix U : 

e{U)^[C DC ... D"-^C], e(L/)GOT„x», (3.33a) 



'B{U)^[B* D*B* ... {D*y'-^B*], G OT„xn, (3.33b) 

and 

§([/) = [C L»C ... D"-^C,B*D*B* ... {D*)"-^B*], 

HU)&mnx2n. (3.33c) 

If the unitary colligation associated with the matrix U is controllable, observable 
or simple, this means that the matrix p.33ap . (j3.33bp or (|3.33cp is, respectively, 
of rank n. 

REMARK 3.4. // one of the matrices (j3.33p has rank n, then its columns ( con- 
sidered as vectors in C" = STtnxij generate the whole space. The columns of these 
matrices are of the form D^C or [D*)^B* , where k takes values in the interval 
[0, . . . , (n — 1)]. It is possible to consider matrices of this kind for k over a larger 
interval. Extending the interval [0, . . . , (n — 1)] does not, however, lead to an in- 
crease in rank for these matrices: The Cayley-Hamilton Theorem tells us that the 
column-vectors, D^C and [D*)^B* with p > n, are, respectively, linear combina- 
tions of the column-vectors D'^C and [D*)^B* with fc G [0, . . . , (n — 1)]. 

DEFINITION 3.14. We say that a unitary matrix U G 2Jt(„+i)x(n+i); 
expressed using the block- decomposition in (I3.28[) . is controllable if 
rank C(C/) = n, observable if rank 23(C/) = n and simple if rank S(?7) = n. 
If the matrix U is both controllable and observable, we say that it is minimal. 

(We note that any one of the matrices p.33p is of rank n if and only if the 
other two have rank n. See Theorem 13.91 ) 

The results of this section on the state space representation of scalar (i.e. 
complex- valued) rational inner functions can be summarized in the following way: 

THEOREM 3.10. ( Rational Inner Functions Equivalence Classes of Unitary Matrices) 

1. Let S{z) be an inner rational function of degree n. Then S{z) can be repre- 
sented in the form: 

S{z) ^ A + zB{Ir,- zD)-^C, (3.34) 
where A, B, C, D are blocks of some unitary minimal matrix U, 

u G9Ji(„+i)x(„+i), mB- 
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2. Let U G $H(n+i)x(n+i) be a unitary matrix with block-decomposition (j3.28[) 
and let the function S{z) be defined in terms ofU by (|3.34[) . Then the function 
S{z) is a rational inner function with degS < n. If the matrix U is minimal, 
then deg S = n. 

3. Let U' G ^{n+i)x{n+i) f^f^d U" G ^(n+i)y.(n+i) ^c Unitary matrices with 
block- decomposition (|3.3ip and let S'{z) and S" {z) be the functions defined 
in terms of U' and U" by: 



S'{z) = A' + zB\In - zD'y'C", S"{z) ^ A" + zS'(/„ - zD")-^C", (3.35) 



// the matrices U' and U" are equivalent, then S'{z) = S"{z). 

If S'{z) = S"{z) and the matrices U' and U" are minimal, then U' and 
U" are equivalent. 

The substance of this theorem can be summarized as follows: 

o There exists a one-to-one correspondence between the set of all rational inner 

functions of degree < n and the set of all equivalence classes of unitary 

matrices in 97l(„+i)x(„+i)- 
o This correspondence can be expressed as a mapping from the set of all rational 

inner functions of degree n onto the set of all equivalence classes of minimal 

unitary matrices in 9Jl(„+i)x(ri+i)- 

For a proof of Theorem 13. 10[ see the Appendix at the end of this paper. 



Applying the Schur algorithm to a given rational inner function s(z) of degree 
n produces the sequence Sfe(z), k = 0, 1, ...n, of rational inner functions with 
sq(z) — s{z) and degsfc(z) = n — i. In particular, Sn{z) = s„ is a unitary constant. 
According to what was stated in Section [3l each of the functions Sk{z) admits a 
system representation, 



in terms of the blocks of some minimal unitary matrix Uk G ^{i+n-k)x{i+n-k)- 



We assume that from the very beginning, the given inner rational function s(z) = 
so{z) is determined in terms of its state space representation, so that the matrix Uq 
is given. The goal is to recursively produce the sequence of matrices Uk representing 
the functions Sfc(z), fc = 1, 2, . . . , ti. The matrix Uk+i, representing the function 
Sk+i{z), is thus constructed from the matrix Uk, representing the function Sk{z). 
In other words, the steps (12. 4p of the Schur algorithm must be described in terms 
of the state space representation (I3.36p . 



The Main Objective of This Paper. 



Sk{z) = Ak + zBk{I - zDk) ^Ck, 



(3.36) 



Uk 



'Ak Bk 
Ck Dk 



(3.37) 
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It should be noted that the unitary matrices in the system representation of 
a rational inner function are determined only up to the equivalence 





Bk 




'1 




'Ak 


Bk 




'l 


0" 


pk 


Dk 









Ck 


Dk 










(3.38) 



where Vk is an arbitrary unitary k x k matrix. So we have to find a rule for 
constructing a matrix Uk+i, which belongs to the equivalence class of matrices 
representing the function Sk+i{z), from an arbitrary element Uk of the equivalence 
class of matrices representing the function Sk{z). 

The Schur algorithm in the framework of system representations 
is described in Section [HI 

Historical Remarks. The definition of a characteristic function was developed 
gradually, starting from the pioneering works of M.S.Livshitz. The first definition 
appeared in [Livl] (for operators for which I — T*T and / — TT* have rank one) 
and in [Liv2] (for the case that these operators have finite rank). M.S.Livshitz and 
those working in the same field, subsequently turned their attention to bounded 
operators T for which T — T* is of finite rank or at least of finite trace. For these 
operators T, a characteristic function was defined in an analogous way and by 
means of this function, a wide-reaching theory for these operators developed. In 
particular, triangular models of non-self-adjoint operators were introduced. See 
[Liv3| , |BrLij , [Br| . In the course of the evolution of the concept of characteristic 
functions, it became clear that it was advantageous to consider, not just non- 
self-adjoint operators, but also more general objects: operator nodes (or operator 
colligations) . The notion of an operator colligation was prompted by physical ap- 
plications of the Livshitz theory of non-selfadjoint operators. (See |BrLi| . [Liv9| 
and references there.) 

B. Sz. Nagy and C. Foias used a different approach to characteristic functions 
in 1962. Their work involved harmonic analysis of the unitary dilation of the 
contractive operator T. Moreover, they simultaneously obtained a functional model 
of T depending explicitly and exclusively on the characteristic function of T. See 
[SzNFol especially Chapter VI] and references therein. 

The version of operator colligations, which appears in Definition 13.11 goes 
back to a remark of M.G. Krein to the work |BrSvl] . In [BrSvl] . the notion of 
a contractive operator colligation (node) was defined as the collection of Hilbert 
spaces H,J-, Q and operators 

To: g ^T, F: J=-^H, G : G ^ H, T : U ^ U , (3.39) 

satisfying the conditions 

/ - TT* = FF*, I - T*T = GG*, 

I - ToTg ^ F*F, I - TqTq = G*G, TG = FTq , (3.40) 
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"To* 


G*' 








'Q' 


-F 


T 




n 




n 



The results presented in the paper |BrSvl| were reported on in a seminar of Krcin's 
in Odessa. In the remark to this talk, M.G.Krein noticed that the conditions (|3.39p - 
(|3.40p mean that the block-operator 



(3.41) 



acting in the appropriate orthogonal sums of Hilbert spaces, is a unitary opera- 
tor. Starting from this remark of M.G. Krein's, mathematicians belonging to the 
Odessa school as well as other mathematicians, defined the operator colligation 

input space . , ,. 

mto the dnect 

state space 

If the spaces have scalar products and the block operator is 



as the block operator acting from the direct sum 



output space 
state space 

a unitary operator with respect to this product, then the operator colligation is 
called an unitary colligation. 

It should be mentioned that the paper [BrSvlj has connections to the theory 
of functional models of contractive operators developed in [Sz NFb] . The definition 
()3.4p of the characteristic function of the colligation (13. 2p - (|3.3p agrees with the 
definition of the characteristic function in [BrSvl . 

The notions of controllability and observability (and minimality) in the set- 
ting of State Space Theory were introduced by R.Kalman in [Kallj . The study of 
controllability and observability of composite systems was first dealt with in [Gil] . 
Under other names, the notion of controllability also appears in the Livshitz theory 
of open systems. See the notions of the simple system and of the complementary 
component in section 1.3 of |Liv9j . (See pages 36-37 of the Russian original, or 
pages 27-29 of the English translation.) 

The fact that every rational matrix-function S can be realized as the transfer 
function of some minimal stationary linear system (which here appears as Theorem 
13. 2p . the uniqueness of the state space representation (Theorem p.4p ) and the 
equality dim Ti. = deg S were all established by R. Kalman in a very general setting. 
These results, as well as many other results, can be found in Chapter 10 of KFA . 
See also Chapter 1 of |Fuhj . 

Some algorithms for the system realization of a given rational function were 
proposed by R. Kalman and his collaborators. (See Chapter 10 of the monograph 
[KFA] and references there.) R. Kalman did not consider questions related to the 
realization of contractive or inner matrix- functions: He developed system theory 
over arbitrary fields rather than over the field of complex numbers. 

An excellent (and short!) presentation of the state space approach to the 
problems of minimal realization and factorization of rational functions can be 
found in [Kaaj . 

Realizations of contractive or inner rational matrix-functions (rational and 
more general) were later considered in the framework of the SzNagy-Foias model 
for contractive operators. These and also more general results can be found in 
many publications now. For convenience, we present some basic facts on system 
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realizations for inner rational functions (scalar) in the Appendix to the present 
paper. 

The state space description of the composite system, which is formed by the 
cascade (or Redheffer) coupling of several state space systems, was dealt with in 
[HeBa] in more generality. We make use of these results, but prefer to derive them 
independently of [HeB a] in the form and in the generality which is most suitable 
for our goal. 



4. Coupled Systems and The Schur Transformation : Input-Output 
Mappings. 

To describe the Schur algorithm using system representations, we must first con- 
sider how the linear-fractional Schur transformation 

So + zcj{z) 



Lu{z) — > s(z), s{z) 



1 



zuj{z)sq 



{sq is a complex number, |so| < 1) 



(4.1) 

can be described in terms of the input-output mappings of linear systems. The 
linear-fractional transform (|4.ip is of the form 



s{z) 



Wii{z)ll>{z) + Wl2(z) 



(4.2) 



W2l{z)uj{z) + W22iz) 

This form of a linear-fractional transform is the most familiar to the classical 
analyst. In the theory of unitary operator colligations, the RedhefFeiH form for 
linear-fractional transforms, i.e. 



(4.3) 



S{z) = Sll(z) + Si2{z)lu{z){I - S22{z)ll!{z)) ^S2l(z) 

is often more convenient. Every linear- fractional transformation of the form (|4.2[) 
can be rewritten in the Redheffer form (|4.3p . but not every transformation in 
Redheffer form can be expressed in linear-fractional form. 

The matrix W(z) — 
(under the appropriate normalizatioro) is 



for the transformation (|4.ip and (|4.2p 



w{z) = ii-wr'^' 



z 
zs^ 



So 
1 



(4.4) 



W{z) in (14. 4|) is not an inner matrix but it is a j-inner matrix: 

j ~W*{z)jW{z) >0, zeB, j -W*{t)jW{t) =0, 

n ' 



t e T, 



where j = 







-1 



^Raymond Redheffer (1921-2005) was a US mathematician working at UCLA. 

*The matrix of the Hnear-fractional transform 1 14.211 is determined only up to the proportionality 

Wiz) X{z)Wiz), where A e C \ {0}. 
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Let us express the fractional-hnear transformation (|4.1 



(|4.3p . where the 2 x 2-matrix-function S{z) — 



Sll(z) Sl2{z) 
S2l{z) S22(z) 



in the RedhefFer form 
is: 



S{z) = 



So 



- \so 



2U/2 



(4.5) 



Unhke W{z), (|4.4p . the matrix-function (|4.5p . is an inner function. 

The transformation in the Rcdhcffer form (I4.3|l admits an interpretation in 
System Theory. We discuss this in more generahty than is needed for our consid- 
erations, which are centered on the linear-fractional Schur transformation. 

Suppose that LSDS and LSDS are two linear stationary dynamical sys- 
tems. In this section, we focus on the input-output mapping and do not touch on 
considerations related to state spaces. 

Let S{z) : f ^ — > 5^ be the transfer matrix-function of the system LSDS^. 
Furthermore, let 

ip{z) = S{z)lp{z) 

be the input-output mapping corresponding to the system LSDS , where (p{z) : 
D £^ is the input signal and tp{z) : I]) ^ £^ is the output signal. Suppose 
now that the exterior space £^ of the system LSDS ^ is the orthogonal sum of the 
subspaces £l and £l: 

£^ = £l © £l . (4.6) 

Equation (|4.6p suggests that the input and output signals be decomposed as fol- 
lows: 



ipiz) 



V'(Z) 



And furthermore that the matrix S{z) be decomposed accordingly: 

Sll{z) Si2{z) 



S{z) 



(4.7) 



(4.8) 



_S2l{z) S22{z) 

So that 

Sll(z) Si2{z) 
S2l(z) S22{z) 

The system LSDS ^ can be considered as a linear stationary dynamical system with 
two input channels, corresponding to the input signals <pi(z) and ip2{z), and two 
output channels, corresponding to the output signals ■0i(^) ^-nd 'ip2iz)'- 







i^2{z)_ 





Vliz) 

f2{z) 



(4.9) 





LSDS^ 




* 


• C 





<— C 



Figure 1 
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Let 

t{z) = uj{z)a{z) (4.10) 

be the input-output mapping corresponding to the system LSDS'', where 17(2;) : 
B is the input signal and t(z) : D is the output signal. The sys- 

tem LSDS^^ can be considered as a linear stationary dynamical system with one 
input channel, corresponding to the input signal <t{z), and one output channel, 
corresponding to the output signal t{z): 



LSDS" 



Figure 2 

Suppose now that 

£'2 = f " (4.11) 

This allows us to 'link' the systems LSDS^ and LSDS^^. We connect the out- 
put channel of the system LSDS with the second LSDS ^ input channel and the 
LSDS^^ input channel with the second LSDS' output channel, as shown in Fig- 
ure 3. 





LSDS' 




• 


» 


















1 


LSDS" 


r 




— > 



Figure 3 

The resulting linear stationary dynamical system LSDS has exterior space 
£l, input signal (pi{z) and output signal tpi{z). The output signal f/'il-^) is linearly 
dependent on the input signal ipi{z): 

Mz) = s{z)Mz), (4-12) 
where s{z) is the transfer function for LSDS. 

We call LSDS the Redheffer coupling of the systems LSDS^ and LSDS". 
We now look to express s{z) in terms of S{z) and ijj{z). The above-described 
connection between the systems LSDS ' and LSDS " can be formally expressed by 
means of the constraints 



^2{z) = t{z), tlj2{z) = a{z) . 



(4.13) 



The Schur Algorithm in Terms of System Reahzations 



27 



EHminating (^2(^), '4'2iz), cr{z), t(z) from the system of linear equations (|4.9p . 
(|4.10p and (|4.13p . we obtain the equation (|4.12p . where s{z) has the form 



S(z) = Sll(z) + Sl2{z)uj{z){l - S22{z)uj{z)) ^ S2l{z) . 



(4.14) 



We now turn our attention to the 'energy relation' associated with the linear frac- 
tional transformation (j4.14|) : uj{z) s{z). 
Equation (|4.9p yields, 



ifilipi + ip^ip2 - Ipllpl - 1p2^2 = [^t '^t] (I - S*S) 

Making the substitutions "01 = st^i, ■02 = o" and ip2 — "-^cr, we obtain 



^l{l-s*s)^i^[^l ^*2]{I~S*S) 



'Pi 

V2 



+ o-*(l -Lu*uj)a , 



(4.15) 
(4.16) 



where 

CT = (1 - S22W)"-^S21</?1, ip2 = w(l - S22i:^)~"^S21t/5l ■ 

It follows from equation (|4l5l) that if 7-5*5 > and l-uj*u; > 0, then l-s*s > 0. 
If / — 5*5 = and 1 — u!*lu = 0, then 1 — s* s — 0. In particular, this brings us to: 



THEOREM 4.1. Let S{z) and uj{z) be rational inner matrix- functions. Further- 
more, let s{z) be given by the Redheffer linear-fractional transform (j4.14p . Then 
s{z) is a rational inner matrix-function. 



We note that the linear-fractional transform, in its classical form (|4.2p . is 
related to another kind of coupling. The relevant connection is shown in Figure 3. 



1 


1 1 
















^1 










SDS 


'P2 


LSDS^ 


1— 










—* 





















Figure 3 

LSDS^ has two input channels with input signals >fi{z) and (p2{z)- LSDS^ also 
has two output channels with output signals ipi{z) and ijj2{z) (with frequency 
representation). The system LSDS^^ has one input channel with input signal <y{z) 
and output signal t{z). We connect the LSDS^^ output channel with the first input 
channel of the system LSDS^ as well as the LSDS^^ input channel with the second 
input channel of the system LSDS' (We assume that the systems are compatible 
with respect to these connections, i.e. that the appropriate subspaces coincide.) 
We consider the second output channel of the system LSDS^ as the input channel 
of the new coupled system LSDS and the first output channel of LSDS^ as the 
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W2l(z) W2l(z) 



output channel of LSDS (Shown in Figure 3.) Let W{z) 
transfer matrix for LSDS^ and lu{z) be the transfer matrix for LSDS^^: 

t{z) = uj{z)a{z) . 



be the 







'wii(z) 


Wl2{z)' 










W2l{z) 


W22{z)_ 




.V^2{z)_ 



The Hnk between the systems LSDS^ and LSDS^^, shown in Figure 3, is described 
by the constraints 

cr{z) ^ (P2{z), T{z)^ifi{z). 

In which the input and the output signals of the system LSDS are denoted by Lp{z) 
and tpiz), respectively: 

= -02(z) , Ipiz) = tpl{z) , 

so that: 

ipiz) = s{z)(p{z) , 

where 

s{z) = {wniz)uj{z) + Wi2(z)) ■ {w2iiz)uj{z) + W22iz)) 

Historical Remark. The coupling of input-output systems having four ter- 
minals, considered in this section (See Figures 1-3), is sometimes called cascade 
coupling. This kind of coupling (as well as related mathematical questions) was 
investigated by R. Rcdhcffcr in |Redl| - [Red5| . Because of this, we use the name 
Redheffer coupling. Redheffer did not consider questions related to cascade cou- 
pling of state space linear systems. These questions were later addressed in [HcBa] 
(Without any reference to Redheffer.) 

The results presented in [HeBaj are more general than here needed. We have 
tailored our approach to the theory of Redheffer coupling in the next two sections 
to fit our needs. 



5. The Redheffer Couphng of Unitary CoUigations. 

As rational inner functions, S{z), ^(z) and s{z) admit system representations as 
characteristic functions of the unitary operator colligations with colligation oper- 
ators U^, U^^ and U, respectively. We now turn to the question of how we might 
express U in terms of the operators and U^^. 

Our approach to this problem will be more general than is here called for, 
our goal being to describe the colligations related to Schur transformations. We 
assume that the unitary colligations corresponding to the systems LSDS ^ and 
LSDS are given. We look to obtain the unitary colligation corresponding to the 
system LSDS, the Redheffer coupling of the systems LSDS^ and LSDS^^. The 
system LSDS ^ is not assumed to be related to the Schur transformation. LSDS ^ 
and LSDS can be generic systems. The only condition imposed on these systems 
is that the exterior space of the system LSDS is identified with the subspace 
£l of the exterior space £i belonging to LSDS^ To avoid technical complications 
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we assume that the exterior and state spaces of the systems LSDS' and LSDS" 
are finite-dimcnsioiial. 

To simpUfy the notation, we denote the matrix entries of the colhgation 
operator corresponding to the system LSDS^, as follows 



where 





ail 


ai2 


u' = 


021 


ail 




_Ci 


C2 






Cp 



hi 
hi 
d 



(5.1) 



The matrix entries for the colligation operator C/^\ corresponding to the system 
LSDS", are denoted as follows: 



a (3 
7 5 



(5.2) 



where 



:£:"^£:", ^:£n^^i 



<5 : 



The linear equations describing the dynamics of the system LSDS ^ are 







ail 


ai2 


hi' 










a2i 


022 


62 




f2{z) 


_z-^h{z)_ 




_Ci 


C2 


d_ 




h{z)_ 



(5.3) 



where 



and h{z) 



- r{z) ■ 




'a 


P 




■a{z)- 






n 


5 




l{z)_ 



are, respectively, the input signal, the output signal and the inner state signal 
corresponding to the system LSDS'. 

The linear equations describing the dynamics of the system LSDS are 



(5.4) 



where cr(z), t{z) and l{z) are, respectively, the input signal, output signal and the 
interior state signal corresponding to the system LSDS^^ 
The constraints 

t{z) = ifiiz) , (j{z) = ipliz) (5.5) 
correspond to the Redheffer coupling of the systems LSDS^ and LSDS". 
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We now aim to eliminate the variables ip2{z)^ ip2iz), cr{z), t{z) from the sys- 



terns (15.31) . (|5.4I) . (|5.5p . To this end, we substitute the expressions aa{z) 
and (t{z) for the variables ^P2{z) and ip2{z) into the equation 

V'2(^) = a2i'^i(2:) + 022932(2) + /3/i(2;) . 
With this we can express <y{z) in terms of ipi{z), h{z) and l{z): 

a{z) = 

= (1 - 022 021 <pi{z) + (1 - 022 &2 Hz) + (1 " 022 a)~^ [3l{z) 
Substituting this expressions for a into (|5.3p . (|5.4p . (|5.5p . we obtain 



Mz) 




- A 


Bi 


B2' 




"Piiz) 


z-^h{z) 












h{z) 


z-H{z) 




.C2 


D21 


D21. 




l{z) 



(5.6) 



(5.7) 



where 
A : £l 



Sl, Bi : Ti.^ Sl , B2 '■ 'hO^ ^ £1, 

C. c I > iy I . ty II > c I 

I . Oi — > n , 02 ■ n — ^ ^1 1 

nt I , n_/ I 7~i . nj II T / I 7~i . T / I 

rt — > rt , JJ\2 - it — * it 1 ^121 ■ it 



The matrix 



U 



Bo 



(5.8) 



A Bi 

Ci D12 
C2 D21 -D21. 

can be expressed using the entries of the matrices U \ (|5.ip . and [/ (|5.2p . as 
follows: 



(1-022 a) ^-[021 &2 022/3]- (5.9) 



The operator U is called the RedhefFer product of the operators Ui and 1/2- 

We again turn our attention to the 'energy relation' associated with the opera- 
tors U^, U^^ and U. Suppose that and U^^ are unitary. Let 
(fi £ £l, f2 & Sl^ h & (T e and I e H^^ be arbitrary vectors. If 
V'l e ^1, V'2 G '?2: ^ S T e iS^^^ and m e are defined by the equali- 
ties 





On 




012/3 




012 a 




Cl 


d 


C2/3 


+ 


C2 a 










5 













i^2 




•^2 


k 




h 



then 
and 



\\k\\' = 

||m|p = 



r 

m 



\\<P2\ 

ii^ir 



(5.10) 
(5.11) 
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For arbitrary ipi, h, I and 

(T = (1 - 022 021 Vl + (1 - 0-22 Oi)~^ 62 + (1 - 022 Ce)~^ (3 , (5-12) 

ip2 = a[{l - 022 a)~^ a2i v?! + (1 - 022 a)^^ &2 + (1 - 022 a)^^ f3) + (31 , (5.13) 
it follows that 

V'2 = Cr, T = (^2 , 

and 

\\i;,\f + \\l\f + \\m\f^\\^^\f + \\h\f + \\l\\\ (5.14) 
According to the definition of the operator U, 



(5.15) 









k 


= U 


h 


m 




I 



Since ipi, h, I are arbitrary, equality (|5.14p means that U is unitary. This operator, 
partitioned into blocks according to (|5.8p . is related to the unitary colligation 
{£,n, U), where £ = £^ , H = T-0- ® H". 

DEFINITION 5.1. The colligation {£,n,U) is called the Redheffer coupling of the 
colligations {£\n\U^) and {£^\ W",?/"). 



THEOREM 5.1. Let S{z) = 



, Lo{z) and s[z) he the characteristic 



Sll(z) Sl2(z) 
S2l(z) «22(z) 

functions of the colligations (5^, 7i^, [/^), {£^^ ,Ti,^\U^^) and their Redheffer cou- 
pling {£,Ti.,U), respectively: 



sii(z) Sl2{z) 
S2i{z) 522(2) 





an 


ai2 




'bi 








+ z 






021 


022 







s{z) = A + z 



{I - zd)-^ [ci C2] , 
[x)(z) — a + zf3{\ — z(5)~"'"7 , 

^ -1 

[Ci C2] 



"Si' 






■ 




'Dn D12 




B2_ 


( 







— z 


D21 D22 


) 



(5.16) 
(5.17) 
(5.18) 



(The notation for the entries of the matrices U\ U^^ and U is taken from (|5.ip . 
(|5.2p and (|5.8p . respectively.) 
Then 



S{z) = Sll(z) + Sl2{z)uj{z){I - S22{z)uj{z)) ^S2l(z) 



(5.19) 
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6. The Inverse Schur Transformation and 
Redheffer Couplings of Colligations. 

We now focus again on the linear-fractional transformation (14. ip in the Redheffer 
form (14. 3p . where uj{z) is a rational inner matrix-function of degree n — 1, so that 
s(z) is a rational inner matrix-function of degree n. 

The function S{z), which appears in (|4.5p is a rational inner function. It is a 
characteristic matrix-function for the operator colligation {£^, Ti^, U^), which we 
now describe. 

The outer space £^ is two-dimensional. We identify £^ with C^. The space 
£^ is considered as the orthogonal sum £^ = £l (B £l, where £l is identified with 
C and £l is identified with C. The orthogonal decomposition £^ — £l® £\ is thus 
the canonical decomposition = C © C. The inner space TL^ is one-dimensional. 
We identify TL^ with C"*^. The colligation operator is defined by the unitary 
3x3= (2 -I- 1) X (2 -I- l)-matrix considered as an operator acting in C'^ = C^: 



with 



So 



B' = 



A' 
"(1 



(6.1) 



l^ol^) 
'So, 



2U/2 



C^=[0 1 ] , D^=[0 



The characteristic function of the colligation (f\7i\C/^) is the matrix-function 
S{z) of the form 



So 



^(l-|soI^) 



2^1/2 



(l-|soP)^/^ 



'ZSq 



(6.2) 



The rational inner function uj{z) of degree n — 1 is the characteristic function 
of the colligation (£",■«",[/"). The outer space £ is one-dimensional and is 
identified with C and the inner space TC^^ is (n — l)-dimensional and is identified 



with C" ^ . The colligation operator U^^ thus acts in C" = C 
the operator U^^ with its matrix in the canonical basis of C 



[/" = 



7 S 



We identify 



(6.3) 



where 



a e OTixi, P e $Hix(„-i), 7 e 2n(«-i)xi, <> e 2n(„_i)x(„-i). 

a is simply a complex number. The matrix U^^ is unitary. The system rep- 
resentation of the function lu{z) is given by: 

uj{z) ^ a + f3{l - zSy'^j . (6.4) 
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In particular, 
The function 



uj{0) = a. 
So + zuj{z) 



1 + zuj{z)so ' 

written as a Redheffer fractional-linear transform, takes the form: 

.s{z) - so + z( 1 - |soP u;{z) ( 1 + zLu{z)J^)-' (1 - \so\')'/' . 



(6.5) 
(6.6) 

(6.7) 



and admits a system realization by means of the operator colligation {£,Ti.,U), 
where {£,H,U) is the Redheffer coupling of the colligations {£^,H^,U^), repre- 
senting the function S{z), and {£^^ , H^^ , U^^) , representing the function uj{z). 

Clearly, £ = C and TC = C". U is the Redheffer coupling of the matrices 
and U^^. Applying formula (|5.9p to and C/^\ we obtain: 

so (l-|so|2)l/2 Oix(„-l)" 

.7(1-^01')^/" -7^ S 
so that U takes the form: 

'A B ' 
C D 



U = 



(6.8) 



(6.9) 



where 



A^So, 5= [(l-|so|2)l/2 Oix(„_i)] 



C = 



D = 


'-a So 


P' 




-ISO 


6 



Clearly, U in §_ 

1 




a(l-|soni/2 
7(1 "kol^)^/^ 
A e aJlixi , B e OTix„ , B e OT„xi, Dem„ 
6.9|) can be expressed as follows: 



(6.10) 



U = 



(ri-l)xl 



Oix(„-i) 
a 13 
7 5 



So 

{l-\so?Y" 

0(„_l)xl 



{i-\s,?Y" 

-So 







(n-l)xl 



Olx(n-l) 
Olx(n-l) 
'-(n-l)x(n-l) 

(6.11) 



Applying Theorem lS.ll to the Redheffer coupling of the colligations (|6.ip . 
and [/", (EH), yields: 

THEOREM 6.1. Let uj(z) be an rational inner matrix-function of degree n — 1 and 
let 



a (3 
7 S 



e OTixi, 13 e 9Jlix(„-i), 7 e fm(„_i)xi, 5 e 9?l(,; 



-l)x(n-l) , 



(6.12) 
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be a unitary matrix so that the system representation (|6.4[) for uj{z) holds. Let sq 
be a complex number with \so\ < 1. Let the function s{z) be defined as the inverse 
Schur transform (|6.6p (using sq and uj(z)) and let the matrix U, 



U 



A B 
C D 



X (ri) 7 



C em^^y,^, D eTl 



(n) X (n) 



(6.13) 



be defined by equation (|6.1ip . 

U is then unitary and yields the system representation of s{z): 

s{z)= A + zB{I - zD)-^C . (6.14) 

Unitary Equivalence Freedom. 

The same function s(z), for which we earlier found a representation using the 
matrix U in (j6.1ip . can also be represented with the help of a matrix having the 
form: 

where V G 97tnxn is a unitary matrix. 

The matrix representing s{z) and which, furthermore, appears as the Redhef- 
fer coupling matrix for the matrices representing S{z) and uj{z), can be considered 
to have fewer 'degrees of freedom' than matrices of the form (|6.15p . The degree 
of freedom for the Redheffer coupling matrix is derived from this same property 
in the Redheffer coupled matrices. The more general form of the matrix, which 
represents the 2 x 2-matrix-function S{z), is the 'transformed' matrix: 



n X n 

V* 



u 



1 

Onxn 



0. 



nxn 

V 



(6.15) 



1 2x2 
1x2 



02x1 

e 



1 2x2 
1x2 



02x1 



I.e. 



U 



A^ 



B^e 



(6.16) 



(6.17) 



where is the matrix from (|6.ip and e is an arbitrary unimodular complex 
number. A more general form of the colligation matrix representing the function 
w{z) is given by: 



U 



11, V _ 



1 

(n-l)xl 



Olx(ri 
V* 



1 

(n-l)xl 



Olx(n-l) 



I.e. 



a 

-1" 



5" 



(6.18) 



(6.19) 



where U (|6.3p , is some nxn unitary colligation matrix representing the function 

uj{z), 

13''= pv, 7"=w*7, S''=v*6v, (6.20) 
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and V is an arbitrary unitary {n — 1) x (ri — l)-matrix. Applying formula (j5.9p to 
the matrices [/^'^ and U^^''", we obtain the Redheffer coupling matrix: 

£(l_|so|2)l/2 Oix(„-l)" 



So 



ae(l-|so|2)i/2 



-a So 



5^ 



(6.21) 



Clearly, 



1 


.0(„-l)xl 



Oix(„_i) 



so e(l-|soP)i/' Oix(„-i) 

e(l-|sop)i/2 -s^ Oix(„-i) 







and finally 



where 



(n-l)xl 

1 ( 

Onxn 







(ti-I)xI 



1 



(n-l)x(n-l). 



v: 



u 



Onx 



Olx(n-l) 



Or 



(6.22) 

(6.23) 
(6.24) 



J(n-l)xl " 

e is an arbitrary unimodular complex number and v is an arbitrary 
unitary (n — 1) x {n — l)-matrix (e and v are the same as in (|6.18p V 

Comparing formulas (|6.15p and (|6.23p - (|6.24p . we see that the matrices C/^' " 
which come from Redheffer coupling of the matrices representing S{z) and oj{z) 
are special. The distinguishing feature of the matrices C/^' " can be summarized as 
follows: 



Among all of the {n + 1) x (n + \)-matrices 



of the form 



so i3n 

(|6.15p . it is precisely those for which the block-matrix entry takes the form 

[£(l-|soP)l/2 Oix(„-l)] , 

where e is an arbitrary unimodular complex number, that can be expressed as 
in (jOSjl - dOl) . 



7. One Step of the Schur Algorithm, Expressed 
in the Language of CoUigations. 

The results from Section[S]can be summarized as follows: Starting from the unitary 

n X ?i-matrix representing a given inner rational matrix-function Ll>{z) of 

7 J 
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so B 
C D 



representing the 



degree n, 

uu{z) = a + zf5{I - zSy^-f , 

we constructed the unitary (n + 1) x (n + l)-matrix 
function s{z): 

s{z) ^ SQ + zB{I - zD)-^C , 

where s{z) is the inverse Schur transform (|6.6p . 

Our goal is not, however, to determine s{z) from uj{z), but instead to start 
with s{z) and determine w(z). We look to describe a step of the Schur algorithm 
when applied to a rational inner function s{z), 



siz) 



1 S{z) - Sq 

uj{z) = -- 



z 1 



So = s(0) , 



s{z)so 
ler word 

representing uj{z), starting from the matrix U representing the 



in terms of system representations. In other words, we would like to find the unitary 

matrix ?1 
7 

function s{z). 



Equation ()6.1ip serves as a heuristic argument. Until now. 



a p 
7 5 



was given 



and U was unknown. Now we assume that the unitary matrix U is given and that 

a [J] , . 

is unknown. We consider (j6.1ip as an equation with respect to 



the matrix 
the matrix 



7 5 
'a 

7 S 



and U as given. Because the second factor on the right-hand 

"a (3' 

7 S 



side of (j6.1ip is a unitary matrix, the solution matrix 



(if it exists) for 



equation (jG.lip is also a unitary matrix. 

For a general unitary matrix U, equation (I6.1ip has no solution 

. \ a P] 
with respect to the matrix 



7 



We know that the block-matrix 



entry B m U = 

|2U/2 



{U as in (|6.1ip ) is necessarily of the form 



So B 
C D 

B= [(l-|so|2)l/2 bix(n-l)] . 

Since the characteristic functions of imitarily equivalent colligations coincide, 
it is enough to find a solution for (|6.1ip with U replaced by some matrix of 
the form (|6l^ : 

2U/2 



1 



(n-l)xl 



Oix(n-l) 

a (3 
7 S 



So 



(1 



- \so\ 

0(ri-l)xl 



1 



Olx(n-l) 
Olx(n-l) 
(n-l)x(n-l). 

(7.1) 
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LEMMA 7.1. Given a unitary (n+ 1) x (n+ \)-matrix U , the unitary n x n-matrix 
V can be chosen such that equation (j7.ip has a solution. 



LEMMA 7.2. Given a unitary {n + 1) x (n + \)-matrix U : 

B^ 



u 



So 
C 



D 



(7.2) 



we can 



find a unitary n x n-matrix Vq such that U^" , given by 



1 



On X n 

takes the form — U'^ , where 

(7° 



On XII 



u 



Onxn 

Vo 



So 
Co 



Bo 
Do. 



(7.3) 
(7.4) 



and the block-matrix entry Bq E SDtixn is 

Bo = [(1- I-SOP)'/' ••• Oix(n-l)] 

PROOF. The row-vectors B and Bq satisfy the condition 

BB*^BoB; (=(l-|soP)) 

The equahty BqBq =1 — |so|^ foUows from the definition of Bq, (|7.4p . The equahty 
so'so + BB* = 1 holds, since the matrix U, ()7.2|) . is unitary. Applying Lemma lSTl to 
the row- vectors B and Bq, we find the unitary nxrt- matrix Vo such that BVo = Bo- 
For every such choice of Vq, the matrix has the form (|7.3p - (l7.4p . 

REMARK 7.1. If n > 1, the matrices and Vq are not uniquely defined. The 
row-vector B of any matrix with V of the form V = ^o[o i where v is an 
arbitrary unitary (n — 1) x (n — l)-matrix is also of the form (|7.4p . 

THEOREM 7.1. Given a unitary (n + 1) x (n + l)-matrix of the form 

So (l-|sol')'/' Oix(n-l)' 



di2 

d22 



ci dii 

C2 ^21 

where sq G C, |so| n > 2, 

CiGOTlixl , diieajlixi , rfi2 e 2tlx(n-l) 

C2 e 9^(n-l)xl , C?21 e 2)t(n-l)xl, '^22 £ 9Jl(n-l) x (n-1) 

the equation 

1 Oix(n-l) 



(7.5) 





(n-l)xl 



a 
7 



13 
5 



So 

il^\so\Y^' 

0(n-l)xl 



-So 

0(n-l)xl 



2U/2 



1 



Olx(n-l) 
Olx(n-l) 

(n-l)x(n-l) 

(7.6) 
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where 

aeajtixi , /3e9Jlix(„_i) , 

7 e 97l(„_l)xl, e 9K(n-l)x(n-l) , 

has a solution with respect to the matrix 

'a P' 



.7 



The solution of this equation can be expressed as 



'a (3' 




-dii So - 


^ci(l- 


|so 


2)1/2 


dl2 


7 6 




-^21 So - 


^ C2 (1 - 




2)1/2 


d22 



(7.7) 



PROOF of THEOREM O We consider equation (|L6j) in further detail. If 
this equation is solvable, then 



So 

Cl 
C2 



(1- 150^)1/2 Oix(„^l)' 

dii di2 

^21 d22 



(1 



|2U/2 







(n-l)xl 



-So 

0(n-l)xl 



Olx(n-l) 
Olx(n-l) 



1 



(n-l)x(n-l). 



1 



(n-l)xl 



Oix(n-l) 

a (3 
7 5 



(7.9) 



Multiplying the matrices on the left-hand side of (|7.9p . we see that their product 

ri on 

is of the form » * * . Since the matrix U , (|7.5p . is unitary, the scalar product of 

L + * 5}i J 

its different rows vanishes. The fact that the first row of this matrix is orthogonal 
to each other row can be expressed as 



Cl 




dii 


C2 


So + 


d2i 



il-\so\y/'^0. 



(7.10) 



The latter equalities mean that the product of the matrices on the left-hand side 
" + * * " 
* * 
* * 

"10 01 

Multiplying out the matrices in 



of (17.911 takes the form 



. Thus, the product of the matrices on the left-hand 



side of 
we obtain 



has the desired form 



* * 
* * 



REMARK 7.2. In view of (|7TTU)) . the solution 
be written as: 



a P 
7 S 



Q.E.D. 

of equation (17. 6p can also 



a (3 




"(1- 


IsoP) 


-l/^Ci d,2 


7 5 




.(1- 


W?) 


-^I^C2 d22 



if |so| < 1, 



(7.11) 
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" a P' 




'-(so) ^dii di2 


7 S 




~{'S^y^d21 d22_ 



if so^O. 



(7.12) 



REMARK 7.3. If n — 1, then there is no room for the matrices (ii2, 1^21, ^22 O'^.d 
13,^,5. In this case U'^, (|7.5p . should be replaced by the matrix: 



so il-\so\Y^' 
ci dii 



where so G C with \so\ < 1, 

CiGOTixi , diiGOTixi, 
and the matrix U^, (|7.7p . should be replaced with: matrix 

= [a] , 

where 

Equation (|7.6p tofces i/ie /orm 



1 

a 



So 

(1-kol^)^/^ 



(i-koH^/^ 



The solution of this equation can be expressed as 

[a] = [-diiSo + ci(l-|so|2)i/2] , 

as well as in the forms: 

H = [(l-|so|')-^/'ci] tf\so\<l, 

and 

[a\^[-{s^)-'dn] if so^O. 



(7.13) 



(7.14) 



(7.15) 

(7.16) 

(7.17) 
(7.18) 



Since both factors on the right-hand side of (|7.6p are unitary matrices, we 
have that is also a unitary matrix. The matrix U'^ in (|7.5p can be considered 
as a matrix of the unitary colligation {£°,TiP, U'^) with outer space £° = C and 
with inner space H'^ — C". The matrix in (|7.7p can, in turn, be considered as 
a matrix of the the unitary colligation {£^,H^, U^) with outer space £^ = C and 
with inner space Ti.^ — C"^^. 

LEMMA 7.3. 

I. If the colligation ,TiP ,U^) is controllable, then the colligation 
(£^,7i^,J7^) is also controllable. 

II. If the colligation {£^ ,T-L^ ,U^) is observable, then the colligation 
(£^,7i^,C/^) is also observable. 
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PROOF. Without loss of generality, we assume that |so| < 1. Otherwise the 
colligation H'^, C/'^) can not be neither controllable nor observable. Our rea- 
soning is based on the equalities 



(7.19) 



'dii 


di2 




'{-So)a 




_d2i 


d22_ 






5 



and 



(7.20) 



Proof of the Statement I. The condition that the colligation (f°,7i°,C/°) be con- 
trollable means that 



V 

0<fc 



'dii 


di2 


k 










d2i 


d22_ 




C2. 


.9^(n-l)xl. 


) 



(7.21) 



And the controllability of the colligation {£^,7i^, U^) can be expressed as: 

V '5'7 = 9K(„_i)xi. (7.22) 



0<fc 



We look to show that (|7?^ follows from (HH} . In view of (TT^ and ([7?^ . we 
can express (I7.2ip as: 



V 

0<fe 



fern 

a 



Let 



"(-so)a 




k 


a 




* 


.(-so)7 


S 




.7. 







OTlxl 

9K(n-l)xl. 



fc = 0, 1, 2 



where fk G 97t(„_i)xi ■ In view of (|7.23p , 

V [fk] =fm(„_i)xi 



0<fe 



Clearly, we have that for every fc = 0, 1, 2, 



/fe-Co,fc'5"7+ ••• +^k-i,kS''-'j + S' 



'1: 



(7.23) 



(7.24) 



(7.25) 



(7.26) 



where S.j.ki < j < fc — 1 , are some complex numbers. Therefore, 



V [^] = V 



0<fe 



0<fe 



We have thus proved Statement I. 
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Proof of the Statement II. The condition that the colhgation be observable 
(f°,7i", C/") can be written as: 

k 



0<k 



dii di2 

^21 d22 



(7.27) 

And the observability of the colligation {£^,Ti,^, U^) can be expressed as: 

V /3'5'^=97Jix(n-i). (7.28) 

0<fc 

We aim to show that (fr28| follows from the formulas (TTf^i . (fTS)) and ((7TT0| . In 
view of (|7.19p . we can express (|7.27p as follows: 

k 



V [1 Olx("-l)] 

0<k 



'(-so)a 



( = [9711x1 97lix(„_i)] ) 



Let 



[l Olx(,i-l)] 



= [* fffe] , A; = 0, 1, 2 



"(-so)a /3" 
.(-so)7 

where gk E 9Jlix(n-i) • In view of (|7.29p . we have that 

V = 3^1x(n-l) • 

0<fc 

Clearly, go = Oix(n-i) and for every fc = 0, 1, 2, . . . 

9k+i = Vo,k f36°+ ■■■ + i^k-i,k PS""-^ + 135'' , 
where < j < fc — 2 , are some complex nmxibers. Therefore, 

V k] = V ■ 



(7.29) 
(7.30) 

(7.31) 

(7.32) 



0<fc 



0<A: 



We have thus proved Statement II. 



Q.E.D. 



The following Lemma is an immediate consequence of Lemma 17.31 
LEMMA 7.4. Let sq G C with \so\ < 1 and [7° be a unitary (n + 1) x (n + 1)- 
matrix of the form (j7.5[) . Suppose that the n x n-matrix , (j7.7p . is related to 
the matrix by equation (j7.6p . Let (S*^ ,TiP ,U'^) and ,Tt^ ,U^) be the above- 
described operator colligations related to the matrices U'^ and U^. If the colligation 
,U^) is minimal, then the colligation ,TC^ ,U^) is also minimal. 

THEOREM 7.2. Let s(z) be a rational inner matrix-function of degree n > 1 (s{z) 
is thus non- constant and \so\ < 1, where sq = s{0)) and let 



, ^ 1 S{Z) - So 

^{z)^--- -r—, so = s(0), 

Z 1-S{Z)SQ 



(7.33) 
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U ^ 



Bo eaJtixn, Co eOT„xi, Do em 



be the Schur transformation of the function s{z). 
Let the unitary matrix U , 

"so Bq 
Co Do_ 

which yields the minimal system representation 

s{z) = So + zBo{I ~ zDo)-^Co , 
have row Bo of the special form 

Bo = [b Oix(„-i)], b>0. 
Then the function lu{z) admits the system representation 
U!{z) = a + zf3{I — z5)~^^ , 

where the unitary n x n-matrix , 

L 7 "J 

aeajlixi, /3 e 97lix(„-i), 7 e 9'J(n-i)xi, (5 e 97t(„_i)x(n-i) 
can be determined from the matrix C/" using 



(7.34) 

(7.35) 
(7.36) 
(7.37) 



a (3 




"(1- 


kol^)' 




7 5 






kol^)- 





(7.38) 



where Cj and djk are the block-matrix entries of the block-matrix decompositions 



Co 



Do = 



dii di2 

d21 d22 



(7.39) 



Cl e ariixi, C2 e 9K(„-i)xi, 

Dii e QJlixl, Di2 e 3Jtlx(n-l): -D2I S 9^(n-l)xlj -^22 G 9K(n-l) X (n-1) 

The unitary colligation associated with the matrix is minimal. 

a (J 



PROOF. The matrix (fZlMl) . the matrix 



7 



(|7.38p . and the number 



So are related by equation (16. lip . According to Theorem 16.11 the function lo{z), 
defined by (|7.37p . and the hmction s{z) are related by the equality (|6.6p . Q.E.D. 



Theorem 17.21 together with Lemma \T72\ describe a step of the Schur algorithm 
in terms of system representations. Before applying the direct Schur transform 
()7.33p . which is a step of the Schur algorithm, we should first 'normalize' the 
colligation matrix U representing the 'initial' function s{z). This normalization 
starts with the matrix U , from which we determine the unitarily equivalent matrix 
, (|7.5p . whose row B^ is of the special form (|7.3p . We then aim to solve the 
equation (|7.6p with respect to the matrix [ " ^ ] • The solution of (|7.6p is given by 

(j7.38p . The unitary matrix [ " ^ ] yields the system representation of the function 
uj{z). 
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It should be emphasized that, in general, the matrix [ " ^ ] is not normalized, 
i.e. its row (3 is not of the form (3 = [* 0ix(„-2)]- To perform the next step of 
the Schur algorithm, we must therefore 'normalize' the matrix [ " obtaining 
the 'normalized' form . We then have to solve the equation of the form (|7.6p . 
where is replaced by C/^, etc. The normalization procedure must therefore 
be performed at every step of the Schur algorithm. This normalization procedure 
is, however, not quite unique. It has some degrees of freedom (See Remark 17.11 ) 
It turns out that we can use these degrees of freedom to make the normalization 
procedure a one-time procedure, so that it might be dealt with during preprocessing 
for the further step-by-step recurrence. In further processing there is then no need 
for normalization and one only has to solve the recurrent chain of equations of the 
form (j7.6p . A one-time normalization of this kind is related to the reduction of 
the 'initial' colligation matrix to the lower Hessenberg form. 

8. Hessenberg Matrices. 
The Householder Algorithm. 

Roughly speaking, the lower (upper) Hessenberg matrix, is a matrix which is al- 
most lower (upper) triangular. The precise definition is: 

DEFINITION 8.1. 

I . We say that a square matrix H is a lower Hessenberg matrix if it has zero-entries 
above the first superdiagonal. If H — | l^ij/c | |o<j,fe<ri; then H is lower Hessenberg 
matrix if hjk = for k>j + l,0<j<n~l. 

II. We say that a lower Hessenberg matrix H ~ \\hjk\\o<j,k<n is special if all 
entries of its first superdiagonal are non-negative: /ij.j+i >0,0<j<n — 1. 

III. We say that a Hessenberg matrix H = ||/ijfc||o<j.fc<n is HL-non-singular if all 
entries of its first superdiagonal are non-zero: 7^0,0<j<n — 1. 

The definition of an upper Hessenberg matrix, special upper Hessenberg ma- 
trix and non-singular upper Hessenberg matrix is similar to Definition 18. II 

DEFINITION 8.2. 

I. We say that a square matrix H is an upper Hessenberg matrix if it has zero- 
entries below the first subdiagonal. If H = ||/ijfe|| o<j,fc<Ti; then H is upper Hessen- 
berg matrix if hjk = for k<j^l,0<j<n — 1. 

II. We say that an upper Hessenberg matrix H = | | |o<j.fe<n is special if all 
entries of its first subdiagonal are non-negative: > , 1 < j < n. 

III. We say that an upper Hessenberg matrix H — ||/ijfc||o<j.A;<n-i is HU-non- 
singular if all entries of its first subdiagonal are non-zero: hjj^i ^ , 1 < j < n. 

Hessenberg matrices were investigated by Karl Hessenberg (1904-1959), a 
German engineer whose dissertation dealt with the computation of eigenvalues 
and eigenvectors of linear operators. 



44 



Bernd Fritzsche, Victor Katsnelson and Bernd Kirstein 



THEOREM 8.1. 

I. Given an (n+ 1) x (n+ \)-matrix M - 
n X n-matrix V such that the matrix , 



1 

Onxl 



Olxn 

V* 



M 



-^•^j,/c||o<i.fe<n7 there exists a unitary 

(8.1) 



1 

Onxl 



Olxn 

V 



is a special lower Hessenberg matrix. 

II. // the matrix M is HL-non- singular, then both matrices and V 
uniquely determined. From the equalities 



1 Oixn 


M 


1 


Olxn 


0„xl V* 


Onxl 


^3 . 



2, 



(8.2) 



where and are special upper Hessenberg matrices, Vi and V2 are unitary 
matrices and the Hessenberg matrix is HL-non- singular, it follows that H2 — 
and Vi^Vx. 

DEFINITION 8.3. Given a square matrix M , a lower Hessenberg matrix H^ to 
which M can be reduced, (j8.ip . is called a lower Hessenberg form of the matrix M. 

THEOREM 8.2. 

I. Given an (n+ 1) x (n+ \)-matrix M — | |-Mj.fe| | o<j,fe<n, there exists a unitary 
n X n-matrix V such that the matrix , 



H"" = 



1 

Onxl 



Olxn 
V* 



M 



1 Olxn 
Onxl V 



(8.3) 



is a special upper Hessenberg matrix. 

II. // the matrix M is HU-non-singular, then both matrices H^ and V are 
uniquely determined. From the equalities 



1 

Onxl 



Olxn 

V* 



M 



1 

Onxl 



Olxi 



1, 2. 



(8.4) 



where H^ and H2 are upper Hessenberg matrices, Vi and V2 are unitary matrices 
and the Hessenberg matrix H^ is HU-non-singular, it follows that H2 = H^ and 

V2=Vi. 

DEFINITION 8.4. Given a square matrix M , an upper Hessenberg matrix to 
which M can be reduced, (j8.3|) . is called an upper Hessenberg form of the matrix 
M. 

THEOREM 8.3. Let U be an {n + 1) x {n+ l)-unitary matrix. 

I. The unitary colligation associated with the matrix U is observable if and only 
if lower Hessenberg form of U is HL-non-singular. 

II. The unitary colligation associated with the matrix U is controllable if and only 
if the upper Hessenberg form of U is HU-non-singular. 



The Schur Algorithm in Terms of System Reahzations 



45 



COROLLARY 8.1. According to Thcorcm l3.9[ the finite-dimensional unitary col- 
ligation is observable if and only if it is controllable. Thus, for a unitary matrix 
U , the lower Hessenberg form of U is H L-nonsingular if and only if the upper 
Hessenberg form of U is HU -nonsingular. 

LEMMA 8.1. Given two row-vectors B' ^[b{b^ . . . . 6/J £ 9Jlix„ and B" = 
[bi b2 ■ ■ ■ ■ 6"] G dJlixn having same norm: 

B'B'* =B"B"* , (8.5) 

there exists a unitary n x n-matrix V such that 

B'V = B". (8.6) 

PROOF of L EM MAO We first consider the question in a more general setting. 
Assume that ^ is a complex Hilbert space with scalar product ( w , w), where {u , v) 
is linear with respect to the argument u and antilinear with respect to v. Let x G i3 
and y G Sj he two vectors such that (x , x) — (y , y) ^ 0. Let ||u|| denote the 

1/2 

norm of the vector u: \\u\ \ = {u , u) ' . Given two vectors x G Sj, y G Sj such that 
Ikll = \\y\\ 7^ 0; our goal is to construct a unitary operator V : S) ^ S) such that 
Vx — y. If the vector y is proportional to the vector x: x — Xy for some A € C, we 
put Vz = Xz Vz e C. This operator is unitary: |A| = 1, because ||a;|| = \\y\\ ^ 0. 
If the vectors x and y are not proportional, we choose A S C, |A| = 1 such that 
A( 2^ , y) > 0. (If ( a; , y) 7^ 0, then this A is unique. If ( x , y} = 0, we can choose 
arbitrary A with |A| = 1.) Let 

Fz = Az-2^fi^^M(Ax-y) VzeSj. (8.7) 
\\x - Ay 11^ 

The vectors 

ei = X + Xy and e2 — x ~ Xy 
are non-zero (x and y are not proportional to one another) and orthogonal: 

(ei,e2) = 0, (8.8) 

because 

(x + Ay , a; - Ay) = (x, x) - AA(y, y) + A(y, x) - A(x, y) 



and (x, x) ~ {y, y), AA = 1, A(x, y) — A(x, y) — A(y, x). From (|8.7p and 
follows that 

Vei = Aei . (8.9) 

From it follows that 

Ve2--~Xe2, (8.10) 

and 

Vz^Xz (8.11) 
Vz £ ^ : {z, ci) — 0, (z, 62) = 0. Therefore the operator V is unitary. Since 

a; =^(61+62), y = -^(61-62) 
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from and ((8?T0)) it follows that Yx = y. 

Let us turn to the proof of the statement of Lemma 18.11 Let be the set of 
all n-row- vectors with complex entries (in other words, = OTixn) and with the 
following scalar product: if u = [ui, ... , m„] and v = [wi, ... , w„] are vectors in 
i^, then their scalar product (m , v) is defined as 

where v* is the Hermitian conjugate of the row-vector v. If H is some n x n-matrix, 
then it generates an operator in ^. This operator maps the row- vector u to the 
row-vector uiJ, where uH is the product of the matrices u and iJ. This operator 
is unitary if and only if H is unitary. 

In the notation of Lemma 18.11 x = B' = [b{ b2 ■ ■ ■ ■ b^^], y — B" = 
[6" 63' • • ■ • ^nl- Thus the matrix V corresponding to the operator (|8.7p takes 
the form 

^^=||^;,fe||i<,,fc<„, (8.12) 

where 

vjk = XSjk - 2(br - Xb^){B' - XB", B' - XB")-\\b'^ - b^) . (8.13) 
and A is such that 

XB'(B")*>Q, |A| = 1. 
5jk is the Kronecker symbol. Q.E.D. 

REMARK 8.1. In the case when the rows B' and B" are real, the matrix V , (|8.12p - 
(|8.13p . is also real. In this case matrices of the form (|8.12p - (|8.13p are known as 
Householder reflection matrices. Householder reflection matrices and the House- 
holder Algorithm (which is based on matrices of this type) are widely used in nu- 
merical linear algebra. See [Wilj . [Strj . [GolV] and |Houj . 

REMARK 8.2. A unitary matrix V satisfying the condition (j8.6p is not unique. 
The process of constructing such matrices (|8.12p - (|8.13p is constructive. 

We will apply Lemma [ST] to the following special situation: Let B' be an 
arbitrary 1 x n-column and B" be of the special form B" = [b" Oix(n-i)]; where 
b" > and thus b" = {B' {B'ff/'^ . For these B',B", the first column of the 
unitary matrix V satisfying (j8.6p is uniquely determined: 

v,,=¥^(B'{Br)-^/\l<]<n. 

The construction of the desired matrix V is thus reduced to the following problem: 
Given the first column of an n x n-matrix, one needs to extend this column to a 
full unitary matrix. The Householder reflection procedure is one way of doing this. 

We use the Householder reflection matrices to reduce an arbitrary matrix to 
a Hessenberg matrix. 
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PROOF of THEOREM [HT] Let M = A/" and let m° ^ be entries of the matrix M°: 



l™i,fcllo<J,fc<n 



(8.14) 



Applying Lemma |8. 11 we choose the unitary matrix Vi G 9Jl,i.„ such that 

["io,i: "io_2, ■ • ■ , "ig jVi = [mj 1, mj_2, ■ • ■ , "Iq J , (8.15) 

where 

Too 1 > 0, mj fc = 0, 2 < fc < n. (8.16) 



(So that mo,i = ([toq ii ; 



'O.li '"0,21 



Vi can be considered as an appropriate Householder rotation, for instance. Let us 
consider the matrix 



1 Oix, 
On X 1 



1 Oixn 

0„xi Vi 



(8.17) 



and let mj j, denote the entries of the matrix M^: 

= l|TO],fcllo<j,fe<n (8.18) 

Clearly, 

mi_o = m°o. (8.19) 

We continue this procedure inductively. We next turn to the inductive step from 
I to I + I. 

Suppose that the matrices G 9Jl„+i_„+i and Vp G dJln-p+i,n-p+i with 
< p < ^ are already known and that the following condition for the entries of 
the matrix Af, 



MP^\\mP\\o<,,k<n. 



are satisfied: 



> 0, m^^ = 0, 3 + 2<k<n, j = 0, 1, 
The matrices Vp, 1 < p < I are unitary and we have 



p-l. 



MP 







(n-p+l)Xf 



px (n— p+1) 
V* 



0(ri— p+1) xp 



0, 



px (n— lp+1) 



Vr 



P 



for every p : p < I . 

We choose the unitary (n — Z) x (n — /)-matrix Vj+i such that 



''"M+lJ ' "'1,1+2^ 



TO; 



: TOL]Vi+i = Ktii, ml+jj: 



,i+2' 



TO 



where 



TOu+i > 0, =0, ; + 2 < fc < n. 



(8.20) 



.21) 



.22) 



.23) 



.24) 



Lemma l8T] ensures that this choice is possible. We then define the matrix Af'+^, 

M'+i = ||TO^.+i||o<,,fc<„ (8.25) 
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as 



M 



i+i 



0(„-;)x(i+l) 







l+lx{n-l) 



Af 



0(n-/)x(i+l) 







(; + l)x(n-i) 



Vi 



+1 



The entries of the matrix Af satisfy the condition 



= 0, j = 0, 1, ... j + 2<fc<n. 



(8.26) 



.27) 



For j — I, condition (|8.24p holds in view of (|8.23p (Ensuring this was our goal in 
choosing the matrix Vi+i as we did.) 

For < j < / — 1, condition (|8.24|) holds, because going from the matrix M' to 
the matrix M'+^ we do not change the rows with indices j : < j < Z — 1: 



0<j</-l, 0< k<n. 



(8.28) 



The equality (|8.28p holds, firstly because the identity matrix of size l + l is the left 



upper corner of the block-matrix 



U! + lx(n-!) 



0(„-!)x(! + l) 



1. / + 1 < A: < n 



and secondly, because 



m^.fe^O Vj,fc: 0<j<l 

(The latter is a consequence of the induction hypothesis (|8.2ip for p = / — 1.) 
The inductive process finishes when we construct the matrix Af„ — A/'^^ for 
l = n-l. 

The matrix V satisfying (|8.ip appears as the product 



V = Vr 
h 



(n-2)x2 



Olx(n- 



(n-2)x2 



02x(n-2) 
1^3 



In-2 
02x(n-2) 



0(n-2)x2 
Vn-1 



5.29) 



According to the above construction, the entries of the matrix H 
= \\hj,k\\o<j,k<n, (ET)), satisfy: 



f^3,k^m^jX^, j<k<n, 



and thus we have: 



(8.30) 
(8.31) 



hj^j+i = > 0, hj^k = 0, j + 2 < k < n, 

Q.E.D. 

The reduction of matrices to the Hessenberg form is a tool often applied in 
numerical linear algebra as a preliminary step for further numerical algorithms. 
See [Wilj . [Strj . [GolV] and other sources in numerical linear algebra. 

The Householder algorithm is implemented in the programming system MATLAB. 
The MATLAB command H=hess(A) reduces the matrix A to the upper Hessenberg 
form H. 

In the next section we discuss the Schur algorithm for rational inner functions 
in terms of the unitary colligation for the system representation of this function. 
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Reducing the coUigation matrix to the upper Hessenberg form is a prehminary step 
for further developing the Schur algorithm in terms of system representations. 

REMARK 8.3. In |KiNe| . the Householder algorithm and the Hessenberg form for 
unitary matrices are used to study the probability measures associated with finite 
Blaschke products via Cayley transform. 



9. The Schur Algorithm in Terms of System Representations. 



We have now finished all necessary preparations and we are well positioned to 
present the Schur algorithm in terms of unitary colligations representing the ap- 
propriate functions. 

Let s(z) be a rational inner matrix-function of degree n > {s{z) is thus 
non-constant and |so| < 1, where so — s(0)) and let 

sq{z) = s{z), Sfc(z), A: = 1, 2, . . . , n , 

be the sequence of rational inner functions constructed according to (|2.4p (deg Sk{z) 
n — fc, so that s„(z) = s„ is a unitary constant.) 
Let 

s{z) ^ A + zB{In- zD)-^C 
be the system representation of s{z), where 

'A B 
C D 



U = 



(9.1) 



(9.2) 



is the matrix of the minimal unitary colligation representing s{z): 
AeMixi, BeMi^n, CeMnxu DeMnxn- {So, A 



So-) 



We first reduce U to the lower Hessenberg form. Let ^ be a unitary n x n-matrix 
such that the matrix (also unitary): 



1 

Onxl 



Olxn 
V* 



U 



1 

Onxl 



Olxn 

V 



J = 1,2. 



is an upper Hessenberg matrix. The block entries of the matrix 

MO ^01 



(9.3) 



(9.4) 



are: 



^" e Mlxl, e OTlxn, C° e fm„xl, e Tlnxn ■ (So, A" = A = so.) 

The unitary colligations associated with the matrices U and U'^ are unitarily equiv- 
alent. The unitary colligation associated with the unitary matrix C/" is therefore 
minimal and represents the function so{z) = s(z): 

so(z) = AO + zBO(/„ ~ zi?")-iC" . (9.5) 

Inductively, we construct the sequence U^, p = 0, 1, . . . , n — 1 of unitary upper 
Hessenberg matrices such that the unitary colligation associated with the matrix 
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jjP = 



(9.6) 



is minimal and represents the function Sp(z), which appears in the p-th step 
of the Schur algorithm. 

For p ~ Q, the representation in (|9.5p holds. We consider the step from p to 
p+1. 

Suppose that W , < p < {n — 1) is a unitary lower iJL- non-singular 
(n — p+1) X (n — p+1) Hessenberg matrix with the block-matrix decomposition: 

'AP BP' 
where 

AP eMi^i, BP e OTix(„-p), CP e M(„„p)xi, dp e 2n(„_p)x(„_p) . 

The unitary colligation associated with the matrix IJP is minimal and represents 
the function Sp{z), which appears in the p-th step of the Schur algorithm: 

Sp{z) =AP + zBP{I„^p - zDP)-^CP . (9.7) 

Let 



CP = 



CI 



DP = 



11 



21 



DP 
DP 



(9.8) 



be the more refined block matrix decomposition of the block-matrix entries CP 
and DP; 

Cf eOTixi, Cf eOT(„_i_p)xi, 

-Dll G OTlixl, -D12 G 37llx(n-l-p), D21 e 9Jl(„_i_p)xl, D22 G 97l(n-l-p)x(n-l-p) • 

Since UP is an upper Hessenberg matrix and also an _ff C/-non-singular matrix, we 
have that BP ^ 0. Because JJp is also unitary, it follows that \Ap\ < 1, i.e. that 

\sp\ < 1 where Sp — Sp(0) . (9.9) 

The row BP is of the form 

s'' = [(l-l=spnl/^Olx(„-p-l)] 

We construct the [n — p) x {n — p)-matrix IJP^^: 

^p+i BP^^' 



UP^ 



(9.10) 



(9.11) 



Af+l GMixl,i?^+^ GOTix( 
where 

(JP+l £)P+^ 



p-l), CP^^ G 9Jl(„_p_i)xl, -D''^^ G 97l(„-p-l)x(n-p-l): 
(l-|,Sp|2)-l/2Cf 



dct 



(1 



spn-v2cp 



-^12 

DP 
^22 



.12) 



To obtain the matrix IJP from IJP^^, one should delete the left column and the 
upper row of the matrix IJP and then recalculate the first column of the resulting 
matrix. The matrix IJP^^ is then an upper Hessenberg matrix. The matrix C/P"*"^ 
is i/L-non-degenerate, because IJP is TJL-non-degenerate and because the first 
superdiagonal of the matrix J7^+^ is a subset of the first superdiagonal of the 
matrix UP. According to Theorem 17.21 (which can be applied to the matrix UP in 
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view of (|9.10p ). the matrix C/^+^ is miitary and the unitary coUigation associated 
with C/^+^ represents the function Sp-|_i(z) appearing in the p + 1-th step of the 
Schur algorithm: 

sp+iiz) = AP+^ + zBP+\ln-p-i - zDP+Y'CP+^ ■ (9.13) 

These considerations do not directly apply when p = ?i — 1. In this case, there is 
no room for B", C", D". However, we can construct 'part' of the matrix (|9.12p : 

^" = (l-|s„_in-i/2c»-i. (9.14) 

(See Remark 17.31 ) The 1 x 1-matrix is unitary, hence it is a unitary constant. 
Clearly, A" = s„, where s„ is the n-th Schur parameter. This completes the 
description of the Schur algorithm for inner rational matrix-functions in terms of 
system representations. Q.E.D. 

REMARK 9.1. It is particularly easy to determine the sequence {DP}p^o,i ... ,n 
of matrices representing the inner operators of the unitary colligations associated 
with the colligation matrices . The matrix makes up the {n ^ p) x (n — p) 
lower-right corner of the matrix D'^ . The inner rational matrix-function s{z) is 
the ratio of two polynomials: 



sp{z) = cp^^^, degxp(^) = n-p, xp(0) = l = 1. (9.15) 



Clearly, 



Xp(z) = det(/„_p - zDf ) , z"-Pxp(l/^) = det (z/„_p - (i?P)*) , (9.16a) 

thus 

Spiz) = Cpdet ((z/„_p - (Z?P)*)(/„_p - zDPy^) . (9.16b) 



10. An Expression for the Colligation Matrix in Terms of the 
Schur Parameters. 

Let s{z) be a rational inner matrix-function of degree n. Let Sp{z), p = 0, 1, . . . , n 
be the sequence of rational inner functions produced by the Schur algorithm from 
the function s{z), as described in (|2.4p . degSp(2:) = n — p. Let , (|9.6p . be 
the colligation matrix of the minimal unitary colligation, which yields the system 
representation (|9.7p of the function Sp. Among all unitary {n — p-\-l) x {n — p-\-l)- 
matrices representing the function Sp we choose a lower Hessenberg matrix U^. 
Such a matrix exists and is unique. 
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The equality (|7.6p . where is taken as the matrix [7° and [7^+^ is taken as 
a p- 

7 6 - 



the matrix [ " ^ ] takes the form 



jjP 



1 Olxl Olx(„-l-p) 

Oixl 

JJP+l 

.0(n_i_p)xl 

(1 - 







(ri — 1— p) X 1 







Olx(n-l-p) 
Olx(n-l-p) 



(ri— 1— p)xl ^(n— 1— p) X (n— 1— p) 

1, 



n-2. 



The latter formula can be rewritten in the equivalent but more convenient form: 



0, 



(n— xp 



px (n— 

UP 



0(n-p)x(p+l) 



0(p+l)x(n-p) 
JJP+l 



02 



xp 



0(n— p— 1) xp 



Opx2 Opx(n-l-p) 



0(n-p-l)x2 



02x(n-p-l) 
^n—l — p 

0<p<n-l. (10.1) 



For p = 0, the matrix on the left-hand side of (110. ip takes the form [J7°] . For 
p = n — 1, the matrix 1/^^-^ takes the form [/" = s„ and the second factor on the 
right-hand side of (jlO.ip takes the form 



In- 







(n-2)x2 



Sn-1 (1- kp„-lP)l/2 



2x(n-2) 



-Sn-1 
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From (fTO?T|) it follows that 



n 



0<p<n-l 



Jpx2 



0, 



px (n—l—p) 



J2xp 











(n—p—l) xp 







(n-p-l)x2 



2x(n-p-l) 



^n—l—p 



(10.2) 

Multiplying the matrices in (|10.2p . we obtain an expression for the entries of the 
matrix which gives us the system representation of the function s{z) in terms 
of the Schur parameters of 5(2:): 

U° = ||'""fc||o<i,fc<n , 



where 



with 



So, 

Aj_iAj_2 
-sj Aj_iAj_2 
A,, 
0, 



•AiAo, l<j<n, 
■ Afc sFT, I < j <n, 



k = Q, 

< j <n - 1 , k = j + 1, 

< j < n - 1, j + 1 < k <n, 

|2n1/2 



(10.3) 



(10.4) 



(10.5) 



^3 — I'^j I 

One can, in the same way, obtain expressions for the matrices of the unitary 
colligations representing the functions Sj, 1 < j < n. 

It should be mentioned that a matrix of the form (|10.3p . (|10.4p appeared 
in the paper [Ger( formula (66')] and was then rediscovered a number times. See 
[Org], [CoS] . [C(m2l Section 2.5], [Siml Chapter 4], pTibt Theorem 2.17]. 



11. On Work Related to System Theoretic Interpretations of the 
Schur Algorithm 

In this section we discuss the connections between the present work and other 
work relating to the Schur algorithm as expressed in terms of system realizations. 
In particular, we discuss the results presented in [AADL] and in [KiNej . 

The paper [AADL] deals with functions of the class S^,, i.e. with the functions 
s{z) meromorphic in the unit disc and possessing the properties: 
1). For every TV and for all points zi, . . . , zn € 1} which are holomorphicity 

points for s, the matrix ||-ft'(2p, 2:<j)||i<p,ij<jV: -^(-ZjC) — i-^(^=)^(C) ^ (Jqcs not 

have more than k negative squares. 
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2). There exists an N and points zi, . . . , zn G ID such that this matrix has 
precisely k negative squares. 

One of the goals of the paper [AADLj is to discuss the Schur algorithm for functions 
from the class S^, in terms of system realizations. In particular, the results of 
|AADL| are applicable to the special cas^l k = 0, in which they can be simplified. 
In our considerations on the algebraic structure of a step of the Schur algorithm 
we will, for the sake of simplicity, restrict ourselves to finite-dimensional systems, 
which correspond to rational inner functions (of, say, degree n) . We now describe 
the relevant result from jAADL| . adopting the notation used there (to make the 
comparison with the results presented in our paper easier). In |AADL| the function 
s{z) is given by 

s{z) ^ so + zB{I - zD)-^C, (11.1) 
where B, C, D are entries of a unitary matrix U , 
\so Si 



U 



C D 



BeOTix„, CeOT„xi, i?eM„x« (11.2) 



It is not explicitly assumed from the very beginning that the entry B of the matrix 
U has the special form (7.36). The matrix U appears as the matrix V, (1.2), in 
[AADLj . Our notation corresponds to that of [AADLj as follows: The objects, 
which appear as 7, v, u, T in formula (1.2) of [AADLj are sq, B* , C, D in our 
formulas (7.34)-(7.35). The state space which is denoted by K. in (1.2) of [AADLj 
is the space Ti — 2Jl„xi (= C") in our paper. 

Let 51(2;) be the Schur transform of the function s{z)^ 

s,{z) = -- '^"^:'' , so = s{0) (11.3) 
z 1 - s{z)so 

(or (7.33) in our paper). According to [AADLj . si{z) is representable in the form 

si{z) = a + zl3{I - zSy^-f , (11.4) 

with 

^ BC, (3= , ^ BDP, 



7= , PC, S^PDP, (11.5) 

where P is the matrix of the orthogonal projector onto the orthogonal complement 
of the vector B* in H, i.e. 

P e Tlnxn, rankP ^ ?i - 1, P"^ = P, P = P* , BP = . 

(Formulas (jll.Sp are the formulas for the entries of the matrix Vi which appear 
on page 11 of [AADLj .) If we would like to represent the image space PTi as the 
space Til = ^{n-i)xi (= C"~^), H = CffiHi, that is, if we would like the matrix 



^ So is the class of contractive functions holomorphic in the unit disc. 
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of the projector P to be of the form P = 
original matrix U with the matrix 





In-l 



then we have to replace the 



U 



1 

Onxr 



0. 



V 



1 Onxr 
Onxri ^ 

where V G 9Jl„xn is a unitary matrix such that 

"0 

In-l 



So Bo 
Co Do, 



(11.6) 



V*PV 

The condition BP = implies the condition Bo 



(11.7) 

= 0. The last equality 



1 

In-l 

means that Bq is of the form Bo = [b Oix{n-i)\ ■ Since the matrix is unitary, 
we have |soP + -Bq-Bq = 1. Therefore, Bo must be of the form 

Bo = [S{1 - \son/^ Oix(„_i)], 

where (5 is a unimodular complex number. The unitary matrix V from pi.6p 
is not unique: In this case, the degrees of freedom are clear, when we consider 

Is 01 



the replacement V ^ V ■ 



V 



where e is an arbitrary unimodular complex 



number and v, v ^{n~i)x{n-i) ^re unitary matrices. Choosing the number e 
appropriately, we can ensure that Bo is of the form 

Bo = [{1 - \son/' Oix(„„i)]. (11.8) 

Let us decompose the matrices Co, Do, which appear as the entries of the matrix 
Uo from PTS]) : 



Co 



Do 



dii di2 

^21 ^22 



(11.9) 



ci e aJlixi, C2 G 5H(„-i)xi, 
Dii e OJlixi, Di2 e 5Hix(„_i), D21 G 2H(,i-i)xi, -D22 e 3^(n-i)x(n-i) • 
The equalities (|11.5p (where B, C, D are replaced by Bo, Co, Do) now take the 



form 



1 



1 



= Ci, /3 = dl2, 



Thus, the matrix 



= C2, 

a ^' 
7 6 



(11.10) 



(11.11) 



from [AADLj . whose entries appear in the representation pi.4p of the function 
51(2;) is the same matrix which appears in our Theorem 7.2 as the matrix (7.38). 
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(The matrix Vi from page 11 of [AADLj can be considered as a coordinate- free ex- 
pression for the coUigation matrix representing the function si{z).) The difference 
between our work and the work |AADLj is not in the results but in the meth- 
ods. The reason for choosing the expression for the coUigation matrix Ui given in 
[AADL] is not fuUy explained. The facts that the matrix C/i is unitary and that 
the matrix Ui represents the Schur transform si of the function s are obtained as 
the result of a long chain of formal calculations. These calculations come across as 
somewhat contrived and do not serve to further our understanding of the subject 
at hand. 

The state system approach is much more transparent. The fact that the 
matrix Ui is unitary is an immediate consequence of our formula (7.6). The fact 
that the matrix Ui represents the function si is a consequence of the interpretation 
of the linear fractional transform (6.6)-(6.7) in terms of the Redheffer coupling of 
the appropriate colligation. 

The paper [KiNe] can also be considered as relevant to our paper. In [KiNej 
the system representation of Schur functions is not considered at all. Neverthe- 
less, in this work the Householder algorithm is used to calculate the sequence of 
numbers, which can be identified with the Schur parameters of the rational inner 
function naturally related to the appropriate unitary matrix. Namely, given a uni- 
tary matrix U G 9K(n+i)x(ji+i)i the measure ji on the unit circle is related to U 
in the following way: fj,{dt) = {E{dt)ei, ei), where E{dt) is the sprectral measure 
of the matrix U and ei = (1, 0, . . . , 0)"^, ei G ^[n+i)y.i- It is assumed that ei is a 
cyclic vector of U . The following equality holds: 

J + zU f l + zt 

-ei = / -^J■idt) (11-12) 



^I-zU Jjl-zt 

The measure /i generates the (finite) sequence of polynomials orthogonal on the 
unit circle. These orthogonal polynomials (<i>fe is monic of degree k) satisfy the 
recurrence relations 

$fe+i(z) = z$fe(z)-sfe$^(^) (11.13) 
n+,{z) ^ z^iz) - SkzMz) (11.14) 

where Sfe, A; = 0, 1, . . . , n are some recurrence coefficients. There are many different 
names for these coefficients. Recently dubbed 'Verblunsky parameters' by Barry 
Simon in [Sim| . On the other hand, the function in (|11.12|) . which we denote by 
p{z) is holomorphic in the unit disc D and has the following properties. 



p(0) = l, piz)+p{z)>0 (zgD). 
Therefore p{z) is representable in the form 

^ \ 1 + 2S(Z) .^x 

P{z)^- (11.15) 

1 — ZS(Z) 

where s{z) is a function holomorphic and contractive in D. Ya. L. Geronimus 
established that the Verblunsky coefficients s{z) in the recurrence relations (|11.13p 



The Schur Algorithm in Terms of System Reahzations 



57 



- (|11.14p are also the Schur parameters of the functions s{z), which appear in 
(jll.lSp . From (|11.12p and (|11.15p it follows that 



J + zU 



zU 



ei = 



1 + zs{z) 
1 — zs(z) 



(11.16) 



In Lemma 3.2 of [KiNej . the following method for finding Schur (=Verblunsky) 
parameters was proposed: First, the given unitary matrix U should be converted 
to Hessenberg form: 



1 

OjlXl 



Olxn 
V* 



u 



1 

Onxl 



Olxn 

V 



(11.17) 



In [KiNej it is claimed that the entries of the (lower Hessenberg) matrix [7° are 
of the form p0.3p - p0.5p . from which the Schur- Verblunsky parameters can be 
found. However, it follows from pi.l6p that 



s{z) ^ A + sB{I - zD)-^C\ 



where 



U = 



A B 
C D 



(11.18) 
(11.19) 



AeMi^i, SeOTixn, C e OT„xi, DeMnxn. 

Thus the formula (|11.18p can be interpreted as the system representation of the 
function s{z). The formula (|11.18p . where s{z) is defined by (|11.16p from U was 
unfamiliar to us, but we do not think that this formula is new. 

In his forthcoming paper [Arl] Yu. M. Arlinskii studied a related question for 
operator-valued Schur functions 8 acting between separable Hilbert spaces. These 
investigations correspond to the operator generalization of the classical Schur al- 
gorithm which is due to Constantinescu (see Section 1.3 in BC .) Yu. M. Arlinskii 
presents a construction of conservative and simple realizations of the Schur al- 
gorithm iterates 0„ of O by means of the conservative and simple realization of 

e. 



Appendix: 

System Realizations of Inner Rational Functions. 

We prove that every complex-valued (i.e. scalar) inner rational function of 
degree n can be represented as the characteristic function of the minimal unitary 
colligation associated with some unitary matrix U € ^{n+i)x{n+i}- Let us denote 
a given rational inner function by S. The operator colligation whose characteristic 
function is S will be constructed as the 'left shift' operator in the appropriate 
space of analytic functions constructed from S. A similar construction appears 
in a paper by B.SzNagy-C.Foias. See [SzNFoi Chapter VI]. The construction of 
B.SzNagy-C.Foias was adapted to unitary colligations in [BrSv2j . 
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1. The space Ks- The most important part of om construction is the Hilbcrt 
space Ks of rational functions. We consider S" as a function defined on the unit 
circle T, i.e. S : T ^ T. As usual, 

i2 ^ {x ; T C, ||x|| < oo}, where ||.t||2 = {x , x) 

and 

{x,y) ^ J x{t) y{t) m{dt) , 

T 

m{dt) is the normalized Lebesgue measure on T. Let and be the Hardy 
subspaces of the space L^: 

Hl = {xeL^ : {x{t) ,t'')=0, fc = -1, ~2, ...}. 
H'i={xeL^ : {x{t) , = 0, fc = 0, 1, 2, . . . } . 

Clearly, 

It is also convenient to consider the functions from and from as functions 
holomorphic in ID) and in D~, respectively. In particular, the evaluation / f(0) 
is defined for every / in and /(oo) = for every / in H^. 
The space Ks is defined as 

Ks = HIq SHl, (A.l) 

where S Hl_ = {S{t) h{t) : h G H\}. Another description of the space Ks is: 

Ks = {x e : X e Hi, xS-^ e Hi }. (A.2) 

It can be shown that the space Ks consists of rational functions whose poles are 
contained in the set of poles of the function S and that dim Ks = deg S. If all zeros 
Zk of S are simple (see then the space Ks is generated by the functions 

{(1 — izfc)~^}i<fc<n. If S has non-simple zeros, the modification of this statement 
is clear. The space Ks is a reproducing kernel Hilbert space. If / G A'5, then 

f{z) = {f{t),Kit,z)), (A.3) 

where the reproducing kernel K(t, z) is: 



2. The left shift operator. The left shift operator T is defined as 

T{m)^{f{t)- }{Q)e{t))-t-^ iorfeHl. (A.5) 

where 

e(t) = l VteT. (A.6) 

This operator is contractive: 

r/f = ii/f -i/(o)p yfeHi. (A.7) 



The Schur Algorithm in Terms of System Reahzations 



59 



The space Ks , considered as a subspace of , is an invariant subspace of the left 
shift operator T. This is evident from the description (jA.2[) of the space Ks- 

3. The construction of the unitary colligation U. The unitary colligation 
{£, H, U) (see Definition 13. 1[) is defined as follows: Let the state space H be the 
space Ks and let the principal operator D be the left shift operator T, (|A.5p . 
restricted to Ks'. 

n^Ks, Df{t) = {f{t)-f{0)e{t))-t-^ y.f(^H. (A.8) 

The equality ||-D/|P + |/(0)P = II/IP, together with the requirement that the 
colligation operator U, (I3.ip - (|3.2p . be unitary, prompts us to define the exterior 
space £ and the channel operator B : TC ^ £ as follows: 

Let f be a one-dimensional Hilbert space which is identified with the vector 
space C over the field C of scalars. We choose the number P — 1 as a basis vector 
in C and will denote this basis vector by 1. Every number e e C, considered as 
an element e of the vector space C, can be presented as e = el, where the factor 
in front of 1 is the same number e, but considered as an element of the field of 
scalars C. 

The channel operator B is: 

{Bfrn ^ fm, yfen. (A.9) 

Equation (|A.7|) ensures that 

\\Bf\\i + \\Dfr^ = \\fr^, yfeH. 

/(O), which appears in (jA.Sp and (|A.9|) . can be represented using the reproducing 
kernel ((X3l) - (fO)) . Let 

k{t) ^1- s{t)^, {^K(t,0)). (A.IO) 

Then 

Bf={f,k)l. (A.ll) 

The operator A : £ ^ £ {as is the case for every operator in £ : dim£ = 1) is of 
the form 

Ae^a{e, 1)1, ee£, 

where a G C. Since the vector 1, which generates £, is orthogonal to H in the 
orthogonal sum £ (BTi-, the unitary property of U implies that 

{Bf,Al) + {Df,Cl)^0 "if en. (A.12) 

Therefore 

a{Bf ,l) + {Df ,C1)^Q VfeH. 

Let us denote 

Cl^l, I en. 

Equation (|A.12|1 means that 

a{f,k) + {Df,l)=0, VfeH. 
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Thus, one should take 

I ^ -a{D*)-'^k, (A.13) 
where D* is the adjoint to the operator D, with respect to the scalar product ( , ). 
We now look to determine the operator D*. The equality 

{Df,g)^{f,D*g) \/f,geH 

means that 

( (fit) - f{0)e{t)) t~\g{t) ) = ( fit), {D*g){t) ) . 
The last equality implies that 

{D*g){t)^P{tg[t)), VgeKs (A.14) 
where P is the orthogonal projector from onto Ks- Clearly, 

{h{t),S{t)) VheKs, 

and 

tg{t) - {tg{t),S{t) )S{t) e Ks V.9 G Kg . 

Therefore, 

P{tg{t)) = tg{t) - ( tg{t), Sit) )S{t) ^g & Kg , 

that is 

{D*g)it) ^tgit) - (tgit),S{t))S{t) VgeH. (A.15) 
From (|A.13p . we obtain 

_ a S{t) - S{0)e{t) 
^ ' ~ S{0) t 
\\Ae\\^ + llCef = 1 gives us \a\ = |S'(0)|. We choose 

a = S{0) 

(Later we see that this is the only possible choice for a.) We set 

lit)='-i^l^^. (A.16) 

(In intermediate steps we assumed that S{0) ^ 0, but this does not appear in the 
final expression (jA.16[) for l{t).) Thus, 

Ae = S{0){e,l)l, Bf = {f,k)l, Ce={e,l)l, 

{Df){t) = {f-{f,k)eit))t-' Vee£,fen. (A.17) 

or 

Ae = eS'(0)l, Bf = f{0)l, (Ce)) = el{t), 

{Df){t)^{f{t)-f{0)e{t))t-\ Ve^eleS, f EH. (A.18) 
From (jA.lSp it follows that the block-operator 



U : 



A B 
C D 



(A.19) 
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is unitary (After the block D was chosen, the other blocks A, B, C were chosen 
to ensure that U he a unitary operator.) The characteristic function Suiz), 

Su{z) = A + zB{I - zD)-^C , (A.20) 

of the colligation U coincides with the original rational inner hmction S{z). This 
can be checked by direct calculation of Sjj{z) using the expression (|A.18[) for blocks 
of the colligation operator U . The expression for the operator (/ — zD)^^, which 
is needed for this calculation, is 

{{I-zD)-'f){t)^ *-^^^\yJ^'\ VfeH. (A.21) 

In what follows we also need the expression for the operator (/ — zD*)~^: 

{{I - zD*r'f){t) ^ ^^^^ ~ ifS-^')iz-')Sit) ^ ^j^^^ (A 

(Since the function fS^^ belongs to H"^ , the evaluation fS^^ — > [f S^^){z^^) is 
defined for z e D.) 

Choosing an orthogonal basis in the (n-dimensional) Hilbert space Ks, we 
realize that the unitary operator [/, (|A.18I) - (|A.19|) . originally constructed as an 
operator acting in a junctional space, is a matrix operator acting in C"+^ = C0C". 

DEFINITION 11.1. The colligation (jATSjl - fATQl) is called the model unitary col- 
ligation constructed from the rational inner function S. 

4. Minimality of the model unitary colligation {£, H, U). We look to prove 
that the model colligation U, (|A.18P - (jA.19|) . is controllable and observable. In view 
of the expression for the channel operator C (one-dimensional), controllability of 
U can be formulated as follows: 

The set of vectors {{I — zD)~^l}z^o generates the space Ks- (A. 23) 

From (|A.16P and (|A.2ip it follows that 

({I - zD)-'l)it) ^ ^W-^(^) . 

t — z 

Let / G be such that 

r Sit) - S{z) — ^ ^ Vz e B (A.24) 
J t- z 

T 

If / e Hi, then / {^m{dt) = Vz G P, hence, / ^^^^ m{dt) = Vz G D. 



The last equahty imphes that f{t)S-^{t) G Hi . (Here we use that S{t) = S-'^{t) 
for t G T.) If als o f {t) S-^{t) G H^ , then f{t)S-\t) = and / = 0. Therefore, 
if the condition (|A.24[) holds for some / G Kg, then / = 0. Controllability of the 
colligation {£, H, U) is thus proved. 
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Observability of this colligation can be proved analogously. According to 
(jA.lip . B* f = {f,e)k. Therefore the observability criterion is reduced to the 
statement: 

The set of vectors {(/ — zD*)~^k}z£o generates the space Kg. (A. 25) 
Using expressions (|A.22[) and (|A.10[) . we obtain: 

Let / £ is such that 

l-S{t)S i(z 1) — VzeB. (A.26) 
1 — tz 

If f{t)S-'^{t) e Hi, then /^^$^ = 0, hence J ^m{dt) = Vz e D and 

T T 

f{t) eHl.lf also f eHl, then / = 0. Therefore if the condition ([XM]) holds for 
some / e Ks, then / = 0. 

5. Uniqueness of simple realization. The uniqueness of the minimal realiza- 
tion is, in fact, a version of a result by M.S. Livshitz, which, in the language of 
M.S.Livshitz, claims that the characteristic function uniquely determines (up to 
unitary equivalence) the operator colligation without complementary component. 

^2 ^2] 
C2 D2 



Let Ui 
into blocks, 



and U2 



be two unitary matrices divided 



where Uj, j — 1, 2, are natural numbers. (A. 27) 

We do not assume that ni — n2- 
Let 

S,iz) ^ A + B^I - zD,)-^C,, i-1,2, 
be the characteristic functions of the unitary colligations associated with the ma- 
trices Ui and U2 respectively. Suppose that 

1. The characteristic functions are equal. 

Si{z) = S2{z). (A.28) 

2. Each of the matrices Ui and U2 is simple in the sense of Definition 13.141 
We prove that under these assumptions the matrices Ui and U2 are equivalent 

in the sense of Definition 13.131 in particular, ni =712- 

To prove this, we have to first of all construct a unitary mapping V of the 
space C"^ onto the space C"^. Assume that the matrices Ui and U2 are simple. 
Let us consider the vectors f^, 5] £ C"^ (= 9Jt„^. xi), j = l,2,0<fc,^: 

= Dp,, g\ = {D*;fB*, j^l,2,0<k,l. (A.29) 
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By the assumption, for each j = 1, 2, the vectors fj, < fc < max(ni, 712) 
generate the space C"^ . The equahty (|A.28|) imphes 

Ai=A2, (A.30) 

and the equahties 

BiDPCi = B2DIC2, < p, (A.31) 

or 

BiD[d'ICi ^ B2DId!^C2, 0<k,l. 
The latter equalities can be interpreted as 

(/i',5l)c"i =(/2',52)c"2, Vfc,;:0<fc, 0<L (A.32a) 
Moreover, the equalities (|A.28[) imply that 

1 - SUOSiiz) ^ 1 - 52*(C)52(z), 1 - Si{z)Sm ^ 1 - S2iz)Sm ■ 
In view of (|3.22p . (|3.23p . the latter equalities imply that 

ci{DiyD{Ci = c;{D;yDPC2, and BiDi{Di)pci ^ B2DI{d;)pc;, 

0<p,q. 

This can, in turn, be interpreted as 

(/f, /i)c"i = (/f, /|)c"2 , and (gf, g1)c"i = (ffa- 9^)c"2 , 

Vp,g : <p, < g. (A.32b) 

From (|A.32p it follows that for arbitrary ak,f3i (such that only finitely many of 
them differ from zero), 

||5]afc/f+^/3,5i||^„, =||^afc/2^- + ^Affyc"2- (A.33) 

Let us define the operator V : C"^ C^^ first as 

Vf^ - /^ Vg'2 ^gl V fc > 0, / > 0, (A.34a) 

and then extend this operator by linearity to all vector columns h G C"^ repre- 
sentable as a finite linear combination of the form h — ^ ctkf2 + S A.92- Thus, 

V{ akf^ + A.92) = E "fc/i' + E 1 • (A.34b) 
If some h G C"^ admits two different representations, say 

/i = E"fe/2 +Ea'5^' /i = E"fc/2 +Ea"5L 

then Vh also admits two different representations: 

F/i = E"fc/' + EA'5i, and F/i = E"fe/i' + E/5"5i- 

However, since X^Q^fc/l + X]A.92 = 0, where = a'^ — a'j,, /?; = /3" — the 
equahty (jA.33p implies that J2 Q^fc/i + J2 Pi9i = 0' i-'^- 

E "'^/i' + E A'ffl = E "'^'/f + E Pi A ■ 
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The definition (fOI]) of V is thus non-contradictory. 

The operator V is defined on the hncar hull of all vectors {/|, 52}fc-' ^^^^ iso- 
metrically maps its definition domain onto the linear hull of all vectors {/i , ffijfc,; • 
If both the matrices [/^, are simple, then these linear hulls are the whole spaces 
C"^ and , respectively. In this case rii = 712 ( n) and 

V*V^In, VV* = /„ (A.35) 

We now prove the intertwining relation 





Bi' 




'l 


0' 




'1 


0' 




'A2 B2 


.^1 


^1. 







V 







V 




C2 D2 



which can be rewritten as follows: 

Ci = VC2, B2 = BiV, (A.36) 

and 

VD2 ^ DiV. (A.37) 

The first of the equalities (|A.36P corresponds to the first of the equalities (jA.29|) 
for fc = (See (jA.34aP for fc = 0.) The second of the equalities (jA.36|) relates to 
the second of equality in (|A.29p for Z = (See (|A.35|) .') 

To check the splitting relation (|A.37p . it is enough to check that 

VD2J2 = DiVf!^ for V/c > 0, (A.38a) 

and 

^£'252 = ^1^92 for > 0, (A.38b) 
The equality (|A.38a[) is an obvious consequence of the definitions of the operator 
V and vectors Indeed, D2/2*^ = /^+\ V f%+^ = f\+^ . On the other hand, 
F/fe = /f, Dif \ = f\+^. Therefore, (IA.38afl h olds. 

Our approach to checking the condition (jA.38bp will be different in the two 
cases / — and / > 0. For I — the equality (jA.38bp takes the form VD2B2 — 
DiVB^. ||X34| for ? = means that VB^ = Bl,so we should check that VD2B^ = 
DiBl . Since DjB* = —CjA*, j = 1, 2, the last equality is equivalent to VC2A2 = 
CiAl. The latter equation is a consequence of the first of the equalities (jA.36[) . 
Thus, (|A.38bp holds for / = 0. 

We check condition (jA.38bp for I > 0. Since the matrices Uj are unitary, we 
have that DjD* = I - CjC*. Thus, (|A.38b|) is equivalent to 

v{i - C2c;){D;y-'B; = (/ - CiC:){Diy''Bi . 

This equation is a consequence of the following three equalities: 

ViD;y-'B; = [Diy-'Bl , (A.39a) 

VC2 = Ci , (A.39b) 

and 

C;iD*y-'^B; = C'*iDiy-^Bl . (A.39c) 
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(|A.39a|) holds, because it can be written as Vg2^^ = which is part of the 

definition (|A.34p of the operator V. (jA.39bp has already been checked: This is 
the first of the relations (|A.36p . (|A.39c|) is the same as (|A.31|) for p = ? - 1. The 
condition (|A.38bp has also been checked for / > 0. 

6. Simple realization is minimal. Let U £ Tl(i^n)x{i+n) be a simple unitary 
matrix. The matrix U is then minimal. Indeed, let S{z) be the characteristic 
function of the unitary colligation associated with U. S{z) is a rational inner 
function, deg5 < n. Let (C, Ks, T) be the model colligation constructed from this 
S. We established that the model unitary colligation is minimal (in particular, 
simple) and that its characteristic function is the function 5, from which it was 
constructed. Both colligations (the original colligation and the model colligation) 
have the same characteristic function and both are simple. Hence, these colligations 
are equivalent. Since the model colligation is minimal, the original colligation is 
also minimal. 
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